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Objective  of  the  Proposal 

This  program,  entitled  “Advanced  Technologies  for  Structural  and  Functional  Optical  Coherence 
Tomography,”  is  a  collaborative  effort  between  investigators  at  the  Massachusetts  Institute  of 
Technology,  Department  of  Electrical  Engineering  and  Computer  Science  and  Research 
Laboratory  of  Electronics;  Dr.  David  Boas,  Associate  Professor  in  Radiology,  Harvard  Medical 
School  at  the  Massachusetts  General  Hospital,  Martinos  Center  for  Biomedical  Imaging;  and  Dr. 
Jay  Duker,  Director,  New  England  Eye  Center,  Professor  and  Chair  of  Ophthalmology,  Tufts 
University  School  of  Medicine.  The  objective  of  this  program  is  to  develop  and  apply  advanced 
optical  coherence  tomography  (OCT)  technologies  for  ultrahigh  resolution  and  functional 
imaging  in  biomedical  applications. 

Accomplishments 

OCT  is  an  emerging  diagnostic  imaging  technology  which  was  developed  by  our  research  group 
in  part  under  this  program  during  the  1990s.  OCT  generates  cross-sectional  images  of  tissue 
microstructure  with  micron  scale  (1-15  pm)  resolution  by  measuring  the  echo  time  delay  and 
magnitude  of  backscattered  light.  OCT  can  thus  function  as  a  type  of  “optical  biopsy,”  enabling 
imaging  of  tissue  with  resolution  approaching  conventional  biopsy  and  histopathology,  but  in 
real  time,  without  the  need  to  remove  and  process  specimens.  During  the  past  contract  period, 
we  have  made  significant  progress  in  a  number  of  areas:  development  of  new  laser  technology 
for  ultrahigh  resolution,  high-speed,  and  three-dimensional  imaging;  new  technology  and 
measurement  techniques  for  functional  OCT  imaging  in  the  retina  and  brain;  in  vitro  imaging 
studies  to  establish  the  feasibility  of  using  OCT  and  Optical  Coherence  Microscopy  (OCM)  for 
optical  biopsy;  and  clinical  studies  in  endoscopy  and  ophthalmology. 

These  key  advances  are  summarized  in  sections  entitled: 

1 .  Optical  coherence  tomography  (OCT)  technology 

2.  Ophthalmic  OCT 

3.  Endoscopic  OCT  studies 

4.  Optical  Coherence  Microscopy 

5 .  Functional  brain  imaging  with  OCT 

A  list  of  publications  supported  under  this  contract  is  appended  at  the  end. 
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1.  Optical  Coherence  Tomography  (OCT)  technology 


Optical  coherence  tomography  (OCT)  is  an  emerging  imaging  technology,  developed  by  our 
research  group  and  collaborators  in  1991  [1].  OCT  enables  imaging  of  the  cross  sectional 
structure  in  biological  tissues  and  materials  by  measuring  the  echo  time  delay  of  backscattered  or 
backreflected  light.  OCT  has  applications  is  a  wide  range  of  clinical  specialties  [2].  OCT  enables 
the  visualization  of  tissue  microstructure  in  situ  and  in  real  time  with  resolutions  in  the  1-10  pm 
range.  The  imaging  depth  of  OCT  in  most  tissues  is  2-3  mm,  determined  by  attenuation  from 
optical  scattering.  OCT  functions  as  a  type  of  optical  biopsy,  enabling  visualization  of  tissue 
structure  without  the  need  for  excisional  biopsy  and  histopathology  [3-6]. 

OCT  is  based  on  low  coherence  interferometry  and  axial  image  resolutions  are  determined  by  the 
coherence  length  of  the  light  source  used  for  imaging.  Clinical  OCT  systems  often  use 
superluminescent  diodes  (SLDs)  that  enable  imaging  with  10-15  pm  axial  resolution.  These 
resolutions  are  typically  insufficient  for  identifying  neoplastic  changes  for  cancer  detection,  or 
tissue  morphological  and  structural  features  for  the  visualization  of  other  pathologies.  Advances 
in  solid-state  lasers  and  nonlinear  fiber  light  sources  have  enabled  the  development  of  ultrahigh 
resolution  and  spectroscopic  OCT  techniques  that  promise  to  improve  tissue  differentiation  and 
image  contrast.  Recent  developments  in  new  detection  techniques  based  on  Fourier  domain 
detection  (known  as  spectral  OCT  and  swept  source  OCT)  enable  very  high  speed.  These  new 
techniques  achieve  dramatic  advances  in  imaging  speed  which  enable  three  dimensional 
imaging.  Tissue  microstructure  can  now  be  visualized  and  rendered  using  methods  similar  to  MR 
imaging,  except  with  micron  scale  resolution. 

1.1  Ultrahigh  resolution  OCT 

Conventional  OCT  imaging  studies  performed  to  date  have  used  axial  resolutions  of  10-15  pm. 
Our  group  pioneered  the  development  of  ultrahigh-resolution  OCT  imaging,  achieving  axial 
resolutions  ~3  pm  in  the  human  eye  and  ~1  pm  in  other  applications  [7-9].  The  axial  resolution 
in  OCT  depends  on  the  coherence  length  of  light  and  is  inversely  proportional  to  the  optical 
bandwidth  of  the  light  source,  k0V AX,  where  Ak  is  the  bandwidth  and  k0  is  the  median 
wavelength.  Ultrahigh  resolution  OCT  continues  to  be  an  active  field  of  research  and  is  closely 
related  to  the  development  of  broadband  light  sources  [2,  10].  This  is  particularly  the  case  for  the 
spectral  region  between  1.0  pm  and  1.5  pm.  This  spectral  region  is  of  great  interest  for  OCT 
because  optical  scattering  is  decreased  at  longer  wavelengths,  enabling  deeper  image  penetration 
depth  in  biological  tissue.  In  addition,  broadband  light  sources  are  essential  for  enabling 
spectroscopic  OCT  imaging  of  functional  parameters  such  as  water  content  and  tissue 
oxygenation  [11]. 

Superluminescent  diode  light  sources  are  commonly  used  in  OCT  because  of  their  compact  size 
and  low  cost.  However,  traditional  superluminescent  diode  light  sources  have  limited  bandwidths 
and  axial  image  resolutions  are  typically  ~10  pm.  Femtosecond  lasers  are  ideal  light  sources  for 
ultrahigh  resolution  OCT  because  they  can  generate  the  extremely  broad  bandwidths  necessary 
for  ultrahigh  resolution  imaging.  Previously,  we  demonstrated  OCT  imaging  with  resolutions  of 
1  pm  at  800  mn  and  5.1  pm  at  1300  mn  in  biological  tissue  by  using  solid-state  mode-locked 


-3  - 


lasers  as  well  as  nonlinear  fiber  sources  [8,  12,  13].  We  have  investigated  compact,  portable  light 
sources  for  ultrahigh  resolution  OCT  imaging  at  wavelengths  of  890  mn,  1300  mn,  and  1500  mn 
to  enable  more  widespread  ultrahigh  resolution  and  spectroscopic  OCT  imaging  to  be  performed 
in  clinical  settings.  Recently,  we  demonstrated  OCT  imaging  with  <5  pm  resolution  in  tissue  at 
1300  mn  and  <3  pm  resolution  at  800  mn  using  continuum  generation  in  a  single  photonic 
crystal  fiber  with  a  parabolic  dispersion  profile  and  two  closely  spaced  zero  dispersion 
wavelengths  [14]. 

1.1.1  Low-threshold  Ti:Sapphire  laser  for  Ultrahigh  Resolution  Ophthalmic  Imaging 

Kerr  lens  modelocked  (KLM)  TiiAFCf  lasers  can  generate  extremely  short  pulse  durations  with 
broad  bandwidths  that  are  particularly  useful  in  biomedical  imaging  [15].  A  standard  Kerr  lens 
modelocked  laser  operating  with  a  5W  pump  can  produce  output  powers  of  500  mW  and 
bandwidths  in  excess  of  150  mn  [16].  Unfortunately,  the  high  cost  of  today’s  femtosecond  lasers 
severely  limits  their  widespread  use.  The  cost  of  femtosecond  TkAfCb  lasers  is  strongly 
dependent  on  the  pump  power  requirements.  Diode  pumped  solid-state  lasers  capable  of 
generating  5  W  can  be  prohibitively  expensive,  while  lasers  generating  several  hundred  mW  are 
considerably  more  affordable. 

We  have  developed  a  low-threshold,  low-cost,  compact  laser  technology.  By  replacing  standard 
mirrors  and  prisms  with  3ld  generation  double-chirped  mirrors  (DCMs),  we  developed  a 
broadband,  portable  laser  source  which  could  be  used  in  the  clinical  setting  for  ultrahigh 
resolution  OCT  imaging  with  3  pm  axial  resolution  [17].  Figure  1.1  shows  the  spectrum  and 
point  spread  function  of  the  low-threshold  source.  Figure  1.2  shows  a  comparison  between 
standard  10  pm  resolution  (a)  and  ultrahigh  3  pm  resolution  (b)  OCT  images  of  the  normal 
human  fovea.  We  developed  a  prototype  ultrahigh  resolution  ophthalmic  instrument  which  was 
used  in  the  ophthalmology  clinic  at  the  New  England  Eye  Center  between  2002  and  2004  and 
performed  studies  on  >1000  patients.  These  studies  provided  new  insights  into  retinal 
pathologies,  especially  alterations  in  photoreceptor  structure. 

1.1.2  Broadband  SLD  for  Ultrahigh  Resolution  Ophthalmic  Imaging 

While  superluminescent  diode  sources  have  proven  to  be  inexpensive  and  stable  for  use  in  a 
clinical  setting,  conventional  SLD  sources  typically  have  limited  bandwidths  and  output  powers. 
Broadband  SLD  sources  are  made  possible  by  multiplexing  two  or  more  diodes  with  different 
center  wavelengths.  These  sources  are  more  inexpensive  and  user-friendly  than  Ti: Sapphire 
sources  but  have  lower  output  powers.  For  applications  where  light  exposure  levels  are  limited 
by  safety  considerations,  such  as  retinal  imaging,  broadband  SLD  sources  provide  sufficient 
power  levels  for  OCT  imaging. 

Recently,  our  group  and  others  have  demonstrated  ophthalmic  retinal  imaging  with  broadband 
SLD  light  sources  [18,  19].  These  studies  were  performed  in  collaboration  with  investigators  at 
Superlum  Diodes,  Ltd.  and  used  a  new  prototype  which  consisted  of  two  spectrally  multiplexed, 
independently  driven  SLD  diodes  operating  at  840  mn  and  920  mn,  respectively.  Figure  1.3 
shows  individual  (a)  and  combined  (b)  SLD  spectra,  along  with  the  point  spread  function  on  a 
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Figure  1.1.  Spectrum  and  point  spread  function  of  the  low-threshold  Ti:Sapphire 
laser  developed  by  our  group  for  OCT  imaging  in  the  ophthalmology  clinic. 
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Figure  1.2.  In  vivo  standard  resolution  OCT  from  the  Zeiss  Stratus  instrument 
(top)  and  ultrahigh  resolution  OCT  images  (bottom)  of  a  normal  human  fovea  at 
approximately  the  same  site.  Resolutions  were  -10-15  pm  axial  x  15  pm 
transverse  (top)  and  -3  pm  axial  x  15  pm  transverse  (bottom)  respectively. 


linear  (c)  and  logarithmic  (d)  scale.  An  SLD  diode  based  on  a  SQW  (GaAl)As  heterostructure 
[20]  and  a  second  recently-developed  SLD  diode  based  on  a  SQW  (InGa)As/(AlGa)As 
heterostructure  with  graded-index  waveguides  [21]  were  used.  The  combined  SLD  output 
provides  155  nm  of  bandwidth  at  a  center  wavelength  of  890  nm  with  >4  mW  of  CW  output 
power.  Compared  to  a  solid-state  laser,  this  broadband  SLD  source  is  quite  compact  and  has  a 
total  footprint  of  only  31  x  26  x  15  cm  including  the  power  supply,  or  31  x  15x5  cm  for  the 
SLD  and  multiplexer  package  alone.  The  total  cost  of  this  light  source  is  about  an  order  of 
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magnitude  less  than  that  of  a  typical  commercial  Ti:Sapphire  femtosecond  laser.  An  image  of  the 
normal  human  macula  taken  with  the  Broadlighter  source  is  shown  in  Figure  1.4.  These  results 
were  significant  because  they  enable  wider  spread  access  to  ultrahigh  resolution  OCT  for 
research  and  clinical  applications. 
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Figure  1.3.  (a)  Individual  output  spectra  of  the  two  superluminescent  diodes,  (b) 
Fiber-coupled  multiplexed  spectrum  of  the  broadband  SLD  source,  (c)  Coherence 
point  spread  function  of  the  broadband  SLD  source,  (d)  Logarithmic  demodulated 
coherence  point  spread  function. 


Figure  1.4.  In  vivo  ultrahigh  resolution  OCT  image  of  the  human  retina  taken 
with  a  broadband  SLD  light  source.  The  axial  image  resolution  in  the  retina  was 
about  3.2  pm.  All  the  major  intraretinal  layers  can  be  clearly  seen  in  this  ultrahigh 
resolution  OCT  image. 
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1.1.3  Cr:Forsterite  for  Endoscopic  Imaging 


Epithelial  cancers  of  the  gastrointestinal  tract,  reproductive  tract,  and  the  respiratory  tract 
comprise  the  majority  of  cancers  encountered  in  internal  medicine.  Many  epithelial  cancers  are 
preceded  by  pre-malignant  changes,  such  as  dysplasia.  Conventional  screening  methods  often 
rely  on  the  gross  morphological  characteristics  of  tissues.  Biopsy  and  histopathology  is  the 
standard  for  the  diagnosis  of  dysplasia  or  carcinoma,  but  can  suffer  from  sampling  errors  and  is 
cumbersome  for  screening  and  surveillance  applications.  OCT  can  function  as  a  form  of  “optical 
biopsy”  which  can  perform  microstructural  imaging  of  tissue  morphology  in  situ  without 
excision.  The  aim  of  this  work  was  to  perfonn  imaging  of  architectural  morphology  of  the 
human  upper  gastrointestinal  tract  to  explore  the  feasibility  of  ultrahigh  resolution  OCT  for  the 
identification  of  morphologies  characteristic  of  early  neoplastic  changes. 

Working  in  collaboration  with  Prof.  Franz  Kaertner,  at  M.I.T.,  we  developed  a  compact, 
femtosecond  Cr4+:Forsterite  laser  combined  with  a  nonlinear  fiber  in  order  to  generate  a  spectral 
bandwidth  greater  than  200  mn  at  a  center  wavelength  of  1250  mn.  The  Cr4+:Forsterite  laser  was 
similar  to  our  previous  systems  [22,  23];  however,  it  used  broadband,  double-chirped  mirrors 
[24]  to  compensate  intra-cavity  dispersion  and  achieve  a  compact,  prism-less  laser  design.  A 
compact  Yb  fiber  laser  was  used  as  the  pump  source.  Figure  1.5(a)  shows  the  optical  spectrum 
generated  by  the  Cr:Forsterite  laser  and  nonlinear  fiber.  The  bandwidth  is  210  mn  FWHM.  Due 
to  bandwidth  limitations  in  the  optical  circulator,  shorter  wavelengths  were  attenuated  and  the 
transmitted  spectrum  was  reduced  to  150  mn  bandwidth  (also  shown).  The  measured  axial  point 
spread  function  shown  in  Figure  1.5(b)  has  a  resolution  of  5  pm  in  air,  close  to  the  calculated 
theoretical  value  of  4.6  pm  for  this  transmitted  bandwidth.  This  corresponds  to  an  axial 
resolution  of  ~3.7  pm  in  tissue,  assuming  an  index  of  refraction  of  -1.37,  and  is  a  three-fold 
improvement  in  resolution  when  compared  to  previous  endoscopic  OCT  systems.  The  system 
detection  sensitivity  was  measured  to  be  102  dB  with  10  mW  of  incident  power. 

Previous  research  suggests  that  OCT  can  image  features  of  architectural  morphology  which  are 
relevant  for  the  diagnosis  and  monitoring  of  early  neoplastic  tissue  changes.  We  have  performed 
in  vitro  surveys  imaging  and  histopathologic  correlations  of  normal  as  well  as  abnormal  human 
tissue  pathologies  including  Barrett’s  esophagus,  adenocarcinoma  of  the  esophagus,  ulcerative 
colitis,  colonic  dysplasia,  adenomatous  polyps,  colorectal  carcinoma,  cervical  intraepithelial 
neoplasia,  cervical  carcinoma,  uterine  dysplasia,  and  endometrial  carcinoma  in  order  to  evaluate 
the  ability  of  OCT  to  differentiate  early  neoplastic  changes  and  other  tissue  abnormalities  with  a 
predisposition  to  malignancy.  In  order  to  validate  this  new  technology  for  ultrahigh  resolution 
endoscopic  imaging  we  perfonned  an  in  vivo  imaging  study  in  a  rabbit  animal  model  [25],  as 
shown  in  Figure  1.6.  Endoscopic  OCT  imaging  was  perfonned  on  the  gastrointestinal,  urinary, 
reproductive,  and  respiratory  tracts  and  correlated  with  histology.  We  also  have  performed 
endoscopic  imaging  of  patients  with  Barrett’s  esophagus. 
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Figure  1.5.  Cr4+:Forsterite  laser  for  ultrahigh  resolution  OCT  at  1300  nm.  (a) 
Optical  spectrum  with  210  nm  FWHM  generated  by  the  Cr4+:Forsterite  laser  and 
self  phase  modulation  in  a  nonlinear  fiber.  Due  to  bandwidth  limitations  in  the 
optical  circulator,  the  transmitted  spectrum  was  reduced  to  150  nm  bandwidth  (also 
shown),  (b)  The  measured  axial  point  spread  function  has  a  resolution  of  5  pm  in  air, 
close  to  the  calculated  theoretical  value  of  4.6  pm  for  the  transmitted  bandwidth. 


Figure  1.6.  (a)  In  vivo  OCT  image  and  (b)  histology  of  rabbit  esophagus  and 
trachea  viewed  intraluminally  from  the  esophagus.  Tracheal  hyaline  cartilage  (he) 
between  the  tracheal  mucosa  and  trachealis  muscle  is  well  defined.  The  image 
demonstrates  the  ability  of  the  endoscopic  OCT  system  to  image  deeply  within 
the  tissue.  Trichrome  stain  was  used  to  highlight  cartilage  and  muscle  layers. 
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1.1.4  All-Fiber  Continuous-Wave  Raman  Continuum  Light  Source 

Working  in  collaboration  with  investigators  from  the  Imperial  College  in  London,  we 
demonstrate  a  new,  high  performance  light  source  for  ultrahigh  resolution  OCT  [26].  Bandwidths 
of  140  mn  at  1300  mn  center  wavelength  with  high  output  powers  of  330  mW  were  generated  by 
an  all-fiber  Raman  light  source  based  on  a  continuous-wave  Yb-fiber  laser-pumped 
microstructure  fiber.  The  schematic  of  the  light  source  and  output  bandwidth  are  shown  in  Figure 
1.7.  The  light  source  is  compact,  robust,  turnkey  and  requires  no  optical  alignment.  In  vivo, 
ultrahigh  resolution,  high-speed,  time  domain  OCT  imaging  with  <5  pm  axial  resolution  is 
possible  using  a  high-power,  continuous-wave,  all-fiber  pump  light  source. 

Microstructure  fibers  typically  have  been  pumped  with  femtosecond  lasers  to  provide  the  peak 
powers  necessary  to  initiate  nonlinear  effects  for  continuum  generation  under  conditions  of 
anomalous  or  near-zero  dispersion.  However,  these  light  sources  can  have  spectral  modulation 
and  excess  intensity  noise  [27].  An  alternative  to  using  high  peak  powers  is  to  increase  the 
effective  nonlinear  interaction  length  of  the  Raman  interaction,  which  is  governed  by  optical 
losses  in  the  fiber  and  dispersive  walk-off  between  the  pump  and  continuum  pulses.  Stimulated 
Raman  scattering  has  been  shown  to  be  the  principle  nonlinearity  for  continuum  generation  by 
using  nanosecond-scale  pump  pulses  [28].  Recently,  the  possibility  of  low  peak  power  and  even 
continuous-wave,  multiwatt  Raman  continuum  generation  in  highly  nonlinear  fibers  has  been 
demonstrated  [29].  Continuous-wave  pumping  of  nonlinear  fibers  can  enable  the  development  of 
robust  and  turnkey  continuum  light  sources  which  require  no  optical  alignment,  enabling  high¬ 
speed,  ultrahigh  resolution  OCT  imaging  in  a  wide  range  of  applications  outside  the  laboratory. 

The  axial  resolution  of  <5  pm  is  comparable  to  that  achieved  using  much  more  complex 
femtosecond  solid  state  sources  such  as  the  Cr:Forsterite  laser.  Output  powers  are  comparable  to 
or  better  than  those  available  from  femtosecond  lasers,  and  more  than  one  order  of  magnitude 
higher  than  superluminescent  diode  sources.  Since  this  light  source  is  all-fiber,  it  requires  no 
alignment  and  provides  completely  turnkey  operation.  The  source  is  extremely  compact, 
measuring  only  25  x  25  x  20  cm,  and  can  easily  be  integrated  into  portable  OCT  systems.  Excess 
noise  is  relatively  high,  but  can  be  reduced  using  dual-balanced  detection.  Further  reduction  in 
noise  should  be  possible  by  reducing  parasitic  feedback  effects  in  the  fiber  splices  and  by 
techniques  such  as  active  seeding  of  the  Raman  process.  Broader  output  bandwidths  should  be 
achievable  using  higher  pump  power  Yb-fiber  lasers  and  different  fiber  geometries.  The 
bandwidth  will  ultimately  be  limited  by  the  water  absorption  of  microstructure  fibers  which 
reduces  the  effective  interaction  length  of  the  Raman  scattering.  Rendering  and  segmentation  of 
sweat  ducts  in  the  human  fingerpad  are  shown  in  Figure  1.8. 

The  current  cost  of  the  source  is  still  relatively  high  due  to  the  cost  of  the  microstructure  fiber. 
However,  since  demand  for  microstructure  fiber  is  increasing,  the  price  of  these  fibers  is 
decreasing  and  can  be  expected  to  approach  that  of  other  specialty  fibers.  The  cost  of  continuous- 
wave  Yb-fiber  lasers  has  also  dramatically  decreased  in  recent  years.  The  cost  of  fiber  Raman 
continuum  light  sources  will  therefore  be  significantly  lower  than  that  of  bulk  solid-state 
femtosecond  laser-based  light  sources.  The  high  perfonnance  and  ease  of  use  of  the  fiber  Raman 
continuum  source  promises  to  enable  a  wide  range  of  new  ultrahigh  resolution,  high-speed,  OCT 
imaging  applications. 
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Figure  1.7.  Left:  Schematic  of  Raman  continuum  light  source.  SMF:  Single-mode 
fiber.  Right:  Output  spectrum.  The  output  was  5.5W  total,  with  2.3  W  in  the  spectral 
range  from  1090  to  1370  mn.  The  output  was  filtered  using  a  special  WDM  coupler  to 
remove  the  pump  wavelength  and  shape  the  spectrum. 


Figure  1.8.  (a)  Rendering  showing  en  face  slice  of  human  fingerpad  created  from 
cross-sectional  images.  Spiraling  sweat  ducts  are  seen  along  the  fingerpad  ridges 
(2.25  x  0.75  x  1.8  mm),  (b)  3D  rendering  and  segmentation  of  sweat  ducts.  Duct 
density  as  well  as  individual  ducts  can  be  assessed. 
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1.1.5  Femtosecond  Nd:Glass  Laser  Light  Source 


High  nonlinearity,  air-silica  microstructure  fibers  or  tapered  fibers  can  generate  a  broadband 
continuum  that  spans  the  visible  to  the  near  infrared  using  femtosecond  pulses.  These  fibers  have 
enhanced  nonlinearity  because  of  dispersion  characteristics  of  these  fibers,  which  shift  the  zero 
dispersion  to  shorter  wavelengths,  and  the  small  core  diameters,  which  provide  tight  mode 
confinement.  High  numerical  aperture  fibers  also  have  been  used  with  femtosecond  Ti: Sapphire 
lasers  to  achieve  bandwidths  of  up  to  200  nm  [30].  Continuum  generation  from  a  femtosecond 
Ti:Sapphire  laser  with  air-silica  microstructured  photonic  crystal  fibers  was  demonstrated  to 
achieve  OCT  image  resolutions  of  2.5  pm  in  the  spectral  region  1.2  pm  to  1.5  pm  [13], 
resolutions  of  1.3  pm  in  the  spectral  region  800  nm  to  1400  nm,  [31]  and  resolution  <1  pm  in  the 
spectral  region  of  550  nm  to  950  nm  [32],  While  high  resolution  imaging  was  demonstrated, 
femtosecond  TiiSapphire  lasers  are  generally  expensive  and  bulky. 

In  this  study,  we  demonstrated  ultrahigh  resolution,  real  time  OCT  imaging  using  a 
commercially  available,  compact  femtosecond  Nd:Glass  laser  that  is  spectrally  broadened  in  a 
high  numerical  aperture  single  mode  fiber  [33].  Ultrahigh  resolution,  high  speed  OCT  imaging 
was  demonstrated  at  1  pm  center  wavelength  with  <5  pm  axial  resolution  in  free  space  (or  -4 
pm  in  tissue)  and  93  dB  sensitivity.  Figure  1.9  shows  the  spectrum  (a)  and  point  spread  function 
(b).  This  system  is  compact  and  well  suited  for  in  vivo  ultrahigh  resolution  OCT  imaging  studies. 
The  1  pm  wavelength  is  important  because  it  has  improved  image  penetration  depths  when 
compared  to  800  nm  and  ultrahigh  resolution  imaging  requires  less  bandwidth  than  at  1.3  pm 
wavelengths  [34].  Figure  1.10  shows  ultrahigh  resolution  images  of  the  hamster  cheek  pouch 
with  the  Nd:Glass  laser  as  the  light  source. 

1.1.6  All-Fiber  Femtosecond  Fiber  Laser  Continuum  Generation 

The  wavelength  region  around  1.4- 1.6  pm  is  of  interest  for  the  characterization  of  optical 
devices  or  spectroscopic  OCT  applications  measuring  water  absorption.  We  investigated  the  use 
of  femtosecond  fiber  lasers  for  real-time,  high-resolution  OCT  imaging  in  the  1.4- 1.7  pm 
wavelength  region  [35].  The  femtosecond  fiber  laser  used  a  ring  cavity  with  an  Er  fiber  gain 
medium  and  passive  modelocking  using  nonlinear  polarization  rotation  as  shown  in  Figure  1.11. 
A  low-noise  supercontinuum  with  38  mW  of  power  and  a  180-nm  bandwidth  is  generated  by  use 
of  a  high-power,  stretched-pulse,  passively  mode-locked  fiber  laser  and  highly  nonlinear  fiber. 
The  spectrum  of  the  fiber  laser  and  the  high-nonlinearity  fiber  are  shown  in  Figure  1.12.  Figure 
1.12  also  shows  the  results  of  imaging  human  skin.  In  vivo  imaging  of  human  skin  is 
demonstrated  with  the  supercontinuum  source  and  a  high-speed  OCT  system.  The  measured 
axial  resolution  was  7.4  pm  full-width-at-half-maximum  (FWHM)  in  air,  which  corresponds  to 

5.6  pm  in  tissue.  This  light  source  is  stable,  compact,  and  self-starting  and  should  be  useful  for  a 
wide  range  of  ultrahigh  resolution,  high-speed  OCT  imaging  applications. 
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Figure  1.9.  (a)  Spectrum  of  the  femtosecond  Nd:  Glass  laser  output  broadened  in 
the  high  numerical  aperture  (NA)  single  mode  fiber,  (b)  Axial  point  spread 
function. 


Figure  1.10.  In  vivo  ultrahigh  resolution  OCT  images  at  1  pm  wavelength. 
Images  were  acquired  at  4  frames  per  second  and  are  2.5  mm  x  1.2  mm  with  1 1 
pm  x  4  pm  (transverse  x  axial)  resolution,  (a,  b)  OCT  of  the  hamster  cheek  pouch 
showing  squamous  epithelium,  connective  tissue  and  a  junction  between  two 
vessels,  (c)  Image  of  the  nail  bed  showing  the  stratum  comeum  and  the  junction 
between  the  epidermis  and  the  dennis.  (d)  Image  of  human  volar  finger  pad 
showing  stratum  corneum  and  sweat  ducts. 
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Figure  1.11.  Fiber  laser  and  nonlinear  fiber  for  continuum  generation  in  the  1500 
nm  region.  The  laser  uses  diode  pumped  Er  doped  fiber  as  the  gain  and  is 
passively  modelocking  using  nonlinear  polarization  rotation.  Broadband 
continuum  is  generated  in  a  nonlinear  fiber. 
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Figure  1.12.  (left)  Output  spectrum  from  the  fiber  laser  and  the  supercontinuum 
generated  in  a  high-nonlinearity  fiber.  The  bandwidth  can  be  broadened  from  70 
to  180  nm.  (right)  High-speed  in  vivo  OCT  image  of  human  skin.  The  top  surface 
is  a  transparent  index  matching  gel  (g).  The  stratum  corneum  (sc),  epidermis  (e), 
sweat  duct  (s),  and  dermis  (d)  are  clearly  visible. 

1.2  High  speed  OCT  using  Fourier  domain  detection 

Conventional  OCT  systems  perform  measurements  of  backscattered  or  backreflected  light  by 
using  an  interferometer  with  a  mechanically  scanned  optical  reference  path  [1,  36,  37].  This 
means  that  light  echoes  with  sequentially  different  delays  are  detected  at  different  times  as  this 
reference  path  length  is  scanned.  Hence,  these  systems  are  known  as  “time  domain”  systems. 
Recently,  novel  Fourier  domain  OCT  detection  techniques  have  emerged  which  do  not  require 
mechanical  scanning  and  achieve  very  high  detection  sensitivities  enabling  OCT  imaging  with  a 
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~15  to  50x  increase  in  imaging  speed  over  standard  resolution  OCT  systems  and  ~100x  over 
conventional  ultrahigh  resolution  OCT  systems.  These  Fourier  domain  techniques  measure  the 
echo  time  delay  of  light  by  Fourier  transforming  the  interference  spectrum  of  the  light  signal  [38, 
39].  Different  echo  time  delays  of  light  produce  different  frequencies  of  fringes  in  the 
interference  spectrum.  Fourier  domain  detection  offers  significantly  improved  sensitivity  and 
imaging  speed  compared  to  time  domain  OCT  [40-42].  Fourier  domain  detection  can  be 
performed  in  two  ways:  spectral  OCT  using  broadband  light  source  and  a  spectrometer  with  a 
linescan  camera  [38,  39,  43-45]  or  swept  source  OCT  using  a  rapidly  tunable,  narrow  linewidth 
laser  source  [1,  46-49]. 

Spectral  and  swept  source  Fourier  domain  OCT  are  especially  promising  for  ultrahigh  resolution 
ophthalmic  imaging  because  they  overcome  imaging  speed  limitations  of  time  domain  OCT. 
Therefore,  it  is  possible  to  use  these  techniques  to  acquire  three-dimensional  maps  of  the  macula 
and  optic  disk  [45,  50].  This  also  enables  cross-registration  of  three-dimensional  data  sets  with 
fundus  photographs  of  the  retina,  for  more  accurate  diagnosis  of  disease  and  evaluation  of 
treatment.  In  addition,  Fourier  domain  OCT  has  the  advantage  of  providing  direct  access  to  the 
spectral  fringe  pattern,  enabling  a  wide  range  of  novel  applications.  Fourier  domain  OCT  can  be 
used  for  absorption  measurement  [51],  Doppler  techniques  can  be  used  to  image  blood  flow[52, 
53],  and  the  complex  Fourier  domain  signal  can  be  directly  measured  to  double  the  axial 
measurement  scan  range  [54,  55]. 

1.2.1  Spectral  domain  OCT 

Spectral  domain  detection  techniques  measure  the  echo  time  delay  of  light  by  using  an 
interferometer  with  a  broadband  light  source  and  measuring  the  spectrum  of  the  interferometer 
output  [40,  43,  51,  54,  56],  Backreflections  or  backscattering  from  the  tissue  at  different  delays 
produces  oscillations  or  fringes  in  the  interference  spectrum.  Increasing  delay  differences 
produce  higher  frequency  oscillations.  The  magnitude  and  delay  of  the  light  from  the  tissue  can 
be  measured  by  Fourier  transforming  the  interference  spectrum.  Spectral  domain  detection  has  a 
sensitivity  advantage  proportional  to  the  number  of  resolvable  elements  in  an  axial  scan.  This 
sensitivity  advantage  is  typically  -20-30  dB,  enabling  dramatic  increases  in  imaging  speeds  [40- 
42,  44,  57].  With  shorter  exposure  times,  spectral  domain  OCT  is  less  sensitive  to  subject  motion 
than  conventional  time  domain  OCT. 

We  designed  and  constructed  a  portable,  high-speed,  ultrahigh  resolution  spectral  OCT  system 
which  can  achieve  axial  image  resolution  of -3  pm  for  imaging  in  the  ophthalmology  clinic.  This 
high-speed  UHR-OCT  prototype  can  be  operated  using  either  a  femtosecond  laser  or  a 
broadband  superluminescent  SLD  as  imaging  light  source.  In  our  early  studies,  we  developed  a 
low  threshold  femtosecond  Ti:  Sapphire  laser  which  was  used  in  the  ophthalmology  clinic  for 
three  years.  However,  recent  advances  in  multiplexed  SLD  light  sources  enable  image 
resolutions  of  3.5  pm.  These  broadband  multiplexed  SLDs  provide  an  easy  to  operate  and 
relatively  inexpensive  light  source  for  clinical  UHR-OCT  imaging.  The  combination  of  high 
image  acquisition  speed  and  ultrahigh  OCT  axial  resolution  yields  significant  improvements  in 
the  clinical  utility  of  high-speed  ultrahigh  resolution  OCT  systems  over  currently  available 
prototype  and  commercial  OCT  systems. 
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Figure  1.13.  Schematic  diagram  of  high  speed  ophthalmic  OCT  instrument  using 
Fourier  spectral  domain  detection.  The  system  uses  a  broadband  light  source  with 
an  interferometer  and  a  spectrometer.  The  interference  spectrum  is  detected  and 
Fourier  transformed  to  obtain  axial  scan  information. 


Figure  1.13  shows  a  schematic  of  the  high-speed,  UHR-OCT  ophthalmic  system  using 
spectral/Fourier  domain  detection  [45].  Light  from  a  broadband  light  source  is  split  between  the 
sample  and  reference  arms  of  an  interferometer.  Light  in  the  reference  arm  is  attenuated  and 
reflected  from  a  stationary  mirror  at  a  fixed  delay.  Light  in  the  sample  arm  is  directed  though  two 
galvanometer- actuated  steering  mirrors  and  relay  imaged  through  the  pupil  onto  the  retina  [37]. 
The  galvanometer  actuated  mirrors  can  scan  the  OCT  beam  across  the  retina  in  any  arbitrary 
pattern  in  order  to  perfonn  cross-sectional  imaging.  OCT  image  data  is  acquired  by  scanning  the 
OCT  beam  on  the  retina  under  computer  control.  The  incident  light  power  on  the  eye  is  750  uW, 
the  same  exposure  used  in  commercial  ophthalmic  OCT  systems  and  consistent  with  ANSI 
safety  standards.  The  spectrum  of  the  interferometer  output  is  detected  using  a  spectrometer 
consisting  of  a  collimating  lens,  transmission  grating,  imaging  lens,  and  CCD  line  scan  camera. 
The  interference  spectrum  data  from  the  camera  was  transferred  to  computer  system  memory 
where  it  was  rescaled  from  wavelength  to  frequency  and  Fourier  transformed  to  generate  axial 
measurements  of  the  echo  delay  and  magnitude  of  light  from  the  retina.  An  example  of  clinical 
imaging  using  the  high  speed,  ultrahigh  resolution  OCT  instrument  with  spectral  detection  is 
shown  in  Figure  1.14. 
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Figure  1.14.  Macular  hole  repair:  fundus  photo  (left  top),  cross-sectional  image 
obtained  by  commercial  instrument.  Zeiss  Stratus  OCT  (top  right),  cross-sectional 
image  obtained  by  high  speed,  ultrahigh  resolution  OCT  instrument.  The  high 
quality  image  consists  of  8100  axial  scans  (transverse  pixels)  with  1024  points 
(axial  pixels)  per  scan  and  was  measured  in  0.35  sec.  NFL  nerve  fiber  layer,  IPL 
inner  plexiform  layer,  INL  inner  nuclear  layer,  OPL  outer  plexiform  layer,  ONL  - 
outer  nuclear  layer,  ELM  external  limiting  membrane,  IS/OS  inner/outer 
photoreceptor  junction,  PR  OS  photoreceptor  outer  segments,  RPE  retinal 
pigment  epithelium. 

Figure  1.14  shows  a  comparison  of  standard  resolution  (10  um  resolution  image  from  the 
commercial  Zeiss  StratusOCT  instrument)  and  high  speed,  ultrahigh  resolution  (from  our 
research  prototype  instrument)  images  from  a  65-year-old  woman  two  years  after  the  surgical 
repair  of  a  full-thickness  macular  hole  in  her  left  eye.  Visual  acuity  had  improved  to  20/25  OS. 
Both  the  standard  resolution  StratusOCT  image  (Fig  1.2b)  and  the  high-definition,  8192  axial 
scan,  ultrahigh  resolution  OCT  image  (Fig  1.2c)  depict  an  epiretinal  membrane  (yellow  arrows). 
However,  the  high-definition,  ultrahigh  resolution  OCT  image  enables  better  visualization  of  the 
membrane  thickness  and  separation  from  the  retina  in  the  nasal  portion  of  the  image.  A  very 
small  central  disruption  of  the  photoreceptor  IS/OS  junction  and  foveal  photoreceptor  outer 
segments  is  visualized  in  the  high-definition  ultrahigh  resolution  OCT  image  (white  arrow).  No 
such  disruption  is  apparent  in  the  standard  resolution  OCT  image.  This  abnormality  in  the 
photoreceptor  signal,  along  with  the  epiretinal  membrane,  may  explain  the  patient’s  slight 
decrease  in  visual  acuity.  Small  photoreceptor  disruptions  visualized  by  OCT  following  macular 
hole  repair  have  been  reported  previously,  and  may  be  correlated  with  visual  acuity. 
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1.2.2  Swept  Source  OCT  imaging 


Although  spectral/Fourier  domain  OCT  systems  can  provide  nearly  optimal  performance  in  low- 
scattering  organs  such  as  the  eye,  these  systems  are  not  well  suited  for  imaging  scattering  tissues. 
Imaging  epithelial  tissues  is  of  great  interest  since  the  majority  of  human  cancers  originate  in 
epithelium  that  lines  hollow  organs  such  as  the  colon,  esophagus,  and  breast  ducts.  In  epithelial 
tissue,  the  high  density  of  cellular  organelles  such  as  mitochondria,  nuclei,  and  cellular 
membranes  results  in  a  highly  scattering  medium  that  limits  the  penetration  depth  of  light  at  800 
nm  wavelengths.  Similarly,  non-biological  objects  such  as  art  samples  also  exhibit  high  optical 
scattering  at  short  wavelengths.  Since  spectral/Fourier  domain  OCT  systems  typically  employ 
silicon-based  CCD  detectors,  they  cannot  operate  at  wavelengths  longer  than  -1000  nm.  For 
OCT  imaging  of  scattering  tissues  and  many  non-biological  samples,  however,  operation  at 
wavelengths  of  1000  -  1300  nm  is  desired  to  reduce  scattering  and  achieve  deeper  image 
penetration  depths. 
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Figure  1.15  Schematic  of  how  swept-source  OCT  operates.  An  RF  beat  signal  is 
generated  by  interference  between  a  sample  beam  and  a  reference  beam.  Optical 
distances  within  the  sample  are  encoded  into  the  RF  beat  frequencies. 

Fourier  domain  detection  approach  using  frequency  swept  lasers  provides  a  powerful, 
complementary  approach  for  spectral  /  Fourier  domain  detection  using  spectrometers  and  line 
scan  cameras.  Swept  source  OCT  systems  typically  employ  a  broadband,  frequency-swept  laser 
source  and  InGaAs  photodetectors  to  perform  OCT  imaging  without  the  use  of  a  spectrometer 
[41,  49,  58-64].  Figure  1.15  illustrates  the  concept  of  swept  source  OCT  imaging  using  a 
frequency-swept  laser  and  a  Michelson  interferometer.  The  frequency  swept  source  light  is  split 
by  a  beamsplitter  into  reference  and  signal  beams.  One  beam  travels  a  fixed  distance  and  reflects 
off  of  a  reference  mirror.  The  other  beam  strikes  the  sample  of  interest,  where  it  is  backscattered 
from  structures  at  various  depths.  Each  backscattering  event  produces  an  attenuated,  time- 
delayed  copy  of  the  incident  frequency  sweep,  which  interferes  with  the  reference  beam.  This 
interference  is  measured  as  a  function  of  time  using  a  high  speed  photodetector.  A  computer 
digitizes  the  interference  fringes  from  the  detector  and  performs  a  Fourier  transformation  to 
generate  an  axial  scan  in  the  OCT  image. 
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Figure  1.16.  RF  beat  signal  formation  using  a  frequency-swept  source  for  OCT 
imaging.  The  detector  outputs  only  the  RF  modulation,  which  is  correlated  to  the 
delay  between  the  sample  and  reference  beams. 

As  shown  in  Figure  1.16,  the  resulting  interference  fringes  contain  RF  beat  frequencies  that 
correspond  to  the  depth  of  each  backscattering  signal.  For  a  single  reflection  in  the  sample  arm,  a 
single-frequency  RF  modulation  will  be  created  where  the  modulation  frequency  is  directly 
proportional  to  the  delay  between  the  sample  and  reference  beams.  For  a  real  material  or 
biological  sample,  each  of  the  numerous  backscattering  signals  creates  a  unique  RF  modulation 
frequency.  After  Fourier  transformation  of  the  RF  beat  signal,  a  single  axial  scan  line  of  an  OCT 
image  is  formed. 
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Figure  1.17.  High-efficiency,  dual-balanced,  fiberoptic  Michelson  interferometer 
used  for  swept  source  OCT  imaging. 

In  practice,  a  fiberoptic  Michelson  interferometer  is  used  for  OCT  imaging.  Figure  1.17 
illustrates  a  typical  imaging  setup  specific  to  swept  source  OCT  systems.  A  small  portion  (~5%) 
of  the  laser  output  is  tapped  off  and  routed  to  a  Mach-Zhender  interferometer,  which  is  used  for 
determining  the  instantaneous  frequency  of  the  sweep.  This  is  necessary  for  resampling  the  OCT 
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interference  fringes  onto  a  uniformly-spaced  frequency  grid  prior  to  Fourier  transformation, 
since  swept  laser  sources  do  not  typically  produce  linear  frequency  sweeps  [65,  66].  A  high- 
efficiency,  dual-balanced  Michelson  interferometer  is  used  to  detect  the  OCT  signal.  In  this 
configuration,  intensity  noise  from  the  laser  source  is  cancelled  out  and  the  effective  signal  level 
is  doubled  by  subtracting  two  out-of-phase  interference  fringes.  The  sample  beam  is  scanned 
over  the  target  by  a  pair  of  XY  mirrors.  A  200  MS/s,  2  channel,  14-bit  A/D  card  is  synchronized 
to  the  laser  sweep  and  is  used  to  record  the  calibration  and  OCT  signals.  A  computer  then 
resamples  the  OCT  signals,  carries  out  a  Fourier  transfonn,  and  produces  a  final  image. 

Several  types  of  frequency-swept  lasers  have  been  developed  for  use  in  OCT  imaging  systems. 
Two  early  studies  were  perfonned  by  our  group  with  bulk  optic  designs,  one  using  a 
Cr:Forsterite  laser  operating  at  1250  nrn  with  a  sweep  rate  of  10  Hz  [46],  and  another  using  an 
external  cavity  diode  laser  operating  at  800  mn  with  a  sweep  rate  of  2  kHz  [47].  An  all-fiber  ring 
laser  employing  a  fiber  Fabry-Perot  tunable  filter  has  been  demonstrated  at  a  sweep  rate  of  200 
Hz  [41].  Hybrid  bulk  optic  and  fiber  cavity  designs  using  a  fiber-coupled  semiconductor 
amplifier  as  a  gain  medium  and  a  diffraction  grating  for  wavelength  selection  have  also  been 
developed.  For  these  designs,  frequency  sweep  rates  of  up  to  115  kHz  have  been  achieved  using 
a  polygon  rotating  mirror  to  sweep  the  operating  wavelength  [64]. 

Since  higher  sweep  rates  allow  correspondingly  higher  OCT  imaging  speeds,  it  is  desirable  to 
increase  the  sweep  rate  to  the  point  where  imaging  becomes  limited  by  detector  bandwidth  and 
A/D  rates.  Recently,  our  lab  has  developed  a  new  class  of  frequency-swept  sources  called 
Fourier  Domain  Modelocked  (FDML)  lasers  [67,  68],  These  lasers  synchronize  the  wavelength 
selective  element,  typically  a  fiber  Fabry-Perot  (FFP)  filter,  to  the  optical  roundtrip  time  of  the 
cavity,  allowing  all  wavelengths  in  the  sweep  to  build  up  and  lase  simultaneously  within  the 
cavity.  Traditional  swept  sources,  in  comparison,  require  lasing  to  build  up  for  each  wavelength 
in  the  sweep  in  series,  limiting  the  maximum  achievable  sweep  rate.  We  have  dsmonstrated 
FDML  lasers  at  sweep  rates  of  up  to  370  kHz,  achieving  OCT  imaging  speeds  of  370,000  A- 
scans  per  second.  Sweep  bandwidths  are  typically  100  -  145  mn  (full  width),  giving  OCT  axial 
resolutions  of  8  -  13  um.  At  these  extremely  high  sweep  speeds,  OCT  imaging  is  limited  by  the 
photodetector  bandwidths  (350  MHz)  and  A/D  rates  (2  MS/s).  For  example,  in  order  to  perform 
true,  real-time  imaging  at  380  kHz,  it  is  necessary  to  use  a  digital  oscilloscope  with  a  5  GS/s 
acquisition  rate  instead  of  a  PC  A/D  card  in  order  to  capture  every  sweep. 

Swept  source  OCT  imaging  has  been  perfonned  in  several  applications  where  rapid  acquisition 
at  1300  mn  is  highly  desirable.  Figure  1.18  shows  an  example  of  OCT  imaging  in  highly 
scattering  human  skin  at  42,000  axial  scans  per  second  and  10  frames  per  second  [67].  Excellent 
penetration  into  the  tissue  is  evident  at  depths  >1  mm.  Due  to  the  high  transverse  line  density 
(4096  lines  per  frame)  and  high  sensitivity  (96  dB)  of  the  swept  source  OCT  system,  the  image 
quality  is  very  high.  Rapid  three  dimensional  imaging  of  epithelial  tissue  is  also  possible  with 
FDML-based  swept  source  OCT  systems.  Figure  1.19  shows  two  orthoslices  of  a  3D  dataset 
taken  from  an  in  vitro  sample  of  hamster  cheek  pouch  fixed  in  formalin.  This  512  x  512  x  200 
pixel  data  was  acquired  at  124,000  axial  scans  per  second,  equivalent  to  242  frames  per  second 
or  1.2  volumes  per  second  [67]. 
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Figure  1.18.  1024  x  4096  pixel  In  vivo  image  of  human  skin  acquired  at  42,000 
axial  scans  per  second.  Good  image  quality  and  penetration  depth  are  observed. 


Figure  1.19.  Two  orthoslices  of  a  512  x  512  x  200  pixel  3D  dataset  taken  from  an 
in  vitro  sample  of  formalin  fixed  hamster  cheek  pouch.  3D  data  was  acquired  at 
124,000  axial  scans  per  second  (1.2  volumes  per  second). 

Another  application  for  high  speed,  swept  source  OCT  is  intravascular  imaging.  In  this 
application,  blood  flow  through  the  vessel  needs  to  be  temporarily  occluded  to  allow  light  to 
reach  the  surface  of  the  vessel.  Occlusion  can  be  perfonned  either  by  temporarily  blocking  the 
vessel  with  an  inflatable  balloon,  or  by  flushing  the  vessel  with  a  saline  bolus.  Both  methods 
result  in  relatively  short  maximum  imaging  times,  from  2-5  seconds  with  the  saline  flush  to  35 
seconds  with  balloon  occlusion.  High  speed  imaging  is  therefore  crucial  to  maximize  the  vessel 
length  that  can  be  investigated  during  the  limited  occlusion  time. 
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We  have  transferred  FDML  laser  technology  to  LightLabs  Imaging  and  this  technology  is  being 
developed  for  their  next  generation  intravascular  imaging  instrument.  Figure  1.19  shows  three 
dimensional  images  of  an  ex  vivo  pig  artery  using  an  FDML  laser  operating  at  45,000  axial  scans 
per  second.  The  imaging  catheter  was  pulled  back  through  the  artery  at  10  mm/s  as  it  rotated  at 
80  revolutions  per  second.  The  FDML  laser  was  operated  at  45,000  axial  scans  per  second, 
giving  4500  OCT  lines  per  mm  in  the  artery.  The  entire  5  cm  artery  segment  was  imaged  in  only 
5  seconds.  This  illustrates  the  potential  of  swept  source  OCT  to  become  a  powerful  tool  for 
applications  which  require  rapid,  survey  imaging  or  large  areas. 


Figure  1.20.  3D  OCT  dataset  taken  from  an  ex  vivo  pig  artery.  Data  was 
collected  at  45,000  axial  scans  per  second,  giving  4500  lines  per  mm  and 
covering  the  5  cm  artery  segment  in  5  seconds. 

2.  Ophthalmic  Optical  Coherence  Tomography 

Since  its  development  in  1991,  OCT  has  become  a  standard  diagnostic  technique  in 
ophthalmology.  Our  group  was  responsible  for  the  invention  of  OCT  and  the  initial  development 
of  OCT  for  ophthalmology.  The  first  report  of  OCT  was  published  in  Science  in  1991  and 
demonstrated  in  vitro  imaging  of  the  human  retina  and  atherosclerotic  plaque  [1].  With  improved 
technology  and  the  availability  of  extensive  clinical  data,  OCT  is  rapidly  becoming  a  standard  of 
care  in  ophthalmology. 

The  current  commercial  OCT  device  (StratusOCT)  has  an  axial  resolution  of  10  pm  and  can 
acquire  an  image  in  ~1.3  seconds.  Recent  research  advances  have  enabled  high-speed,  ultrahigh 
resolution  (UHR)  OCT,  achieving  axial  resolutions  as  fine  as  2  pm  and  acquiring  retinal  images 
in  -0.02  seconds  [45,  50,  69,  70].  Our  current  research  focuses  on  developing  high-speed, 
ultrahigh  resolution  OCT  for  clinical  and  fundamental  research  applications. 

In  many  retinal  diseases,  clinical  outcome  is  significantly  improved  by  detecting  disease  during 
its  early  stages.  Early-stage  retinal  disease  is  often  characterized  by  small  changes  in 
architectural  morphology  or  physiology.  There  is  a  need  to  comprehensively  screen  the  retina  for 
such  changes,  and  quantitative,  objective  markers  for  disease  progression  and  response  to 
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treatment  are  required.  Our  high-speed  UHR-OCT  prototype  enables  three-dimensional  imaging 
of  the  retina,  acquiring  comprehensive,  micron-scale  structural  information  necessary  for  early 
detection  of  retinal  disease.  We  have  developed  a  prototype  instrument  which  is  in  use  at  the 
New  England  Eye  Center  for  clinical  studies  imaging  a  range  of  retinal  pathologies.  In  addition 
to  performing  qualitative  cross-sectional  studies  of  retinal  diseases,  our  group  is  developing 
algorithms  for  quantitative  analysis  of  OCT  images  and  the  creation  of  age-matched  normative 
databases.  The  analysis  of  image  information  and  the  extraction  of  quantitative  metrics  is  an 
important  step  toward  developing  clinically  useful  diagnostic  protocols. 

Another  area  of  research  is  structural  and  functional  OCT  imaging  in  the  murine  (rat  and  mouse) 
eye.  The  murine  eye  is  important  due  to  its  structural  similarity  to  the  human  eye  and  relative 
ease  of  genetic  manipulation,  making  it  an  important  surrogate  model  for  a  number  of  retinal 
diseases.  Our  initial  demonstration  and  subsequent  investigation  of  in  vivo  functional  imaging  in 
the  murine  eye  may  pave  the  way  to  eventual  applications  of  functional  OCT  in  the  human 
retina. 

2.1  Ophthalmic  ultrahigh  resolution  OCT  instruments 

The  OCT  instruments  that  we  used  for  ophthalmic  imaging  are  based  on  either  time  domain 
detection  [1]  (Figure  2.1  a)  or  spectral/Fourier  domain  detection  [38,  39]  (Figure  2.1  b). 
Ultrahigh  resolution  imaging  (~2  -  3.5  pm)  was  achieved  with  a  femtosecond  Ti: Sapphire  laser 
light  source  [8,  71,  72]  (bandwidth  of  150  nm  with  central  wavelength  of  820  nm)  or  broadband 
superluminescent  diode  light  source  [73]  (bandwidth  of  90  nm  centered  at  840  nm).  High  speed 
data  acquisition  rates  (up  to  25,000  axial  scans  per  second)  were  attained  using  spcctral/Fouricr 
domain  detection. 


Rt-’fciejice  ami 


Figure  2.1.  Schematics  of  OCT  systems  for  ophthalmic  imaging:  (a)  Time  domain  OCT  system, 
(b)  spectral/Fourier  domain  OCT  system,  (c)  Retinal  scanning  module  (RSM).  FC  fiber  coupler, 
PC  polarization  controller,  CF  collimating  lens,  RSDF  rapid  scanning  delay  line,  RSM  retinal 
scanning  module,  PD  photodiode,  DC  dispersion  compensation,  NDF  neutral  density  filter,  DG 
-  diffraction  grating,  SOF  spectrometer  objective  lens,  CCD  line  scan  CCD  camera. 
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To  provide  the  best  visualization  of  retinal  pathologies,  we  developed  and  investigated  different 
image  acquisition  protocols.  Imaging  protocols  initially  proposed  by  our  group  in  the  mid  1990s 
have  been  implemented  in  commercially  available  devices  [37,  74].  However,  to  take  advantage 
of  the  enhanced  performance  of  high-speed  spectral/Fourier  domain  OCT,  new  imaging 
protocols  are  required  [44].  One  protocol  acquires  a  small  number  of  6  mm  high-definition 
images  (8192  A-scans  x  3  B-scans).  It  is  especially  useful  for  patients  with  opaque  media  or 
other  conditions  that  result  in  a  low  OCT  signal,  as  transverse  pixel  averaging  can  be  used  to 
increase  the  signal-to-noise  ratio.  A  second  protocol  (Figure  2.2  a)  acquires  a  series  of  high 
transverse  pixel  density  images  in  a  6  mm  x  6  mm  raster  pattern  (2048  A-scans  x  21  B-scans).  A 
third  protocol  (Figure  2.2  b)  acquires  three-dimensional  OCT  (3D-OCT)  data.  It  achieves 
comprehensive  retinal  coverage,  taking  measurements  on  a  6  mm  x  6  mm  area  with  a  spacing  of 
12  pm  x  33  pm  between  axial  scans.  The  large  set  of  cross-sectional  images  can  be  useful  for 
tracking  pathologies  in  three  dimensions,  or  detecting  small  focal  pathologies.  A  fourth  protocol 
(Figure  2.2  c)  consists  of  24  radial  scans  with  1500  lines,  centered  at  the  fovea.  This  protocol  can 
be  used  to  map  intraretinal  layers  and  samples  most  densely  in  the  fovea,  the  region  most 
important  for  visual  acuity. 


Figure  2.2.  Examples  of  OCT  scan  protocols:  (a)  series  of  high  transverse  pixel 
density  images  in  a  6  mm  x  6  mm  raster  pattern  (acquisition  time:  1.8  seconds), 

(b)  three-dimensional  OCT  data  (acquisition  time:  3.7  seconds),  (c)  radial  scans 
(acquisition  time:  1 .4  seconds). 

2.2  Studies  on  normal  subjects 

High-speed  UHR-OCT  has  a  number  of  advantages  including  improved  image  quality, 
preservation  of  retinal  topography,  improved  retinal  coverage.  Images  obtained  using  this 
technique  show  new  features  at  the  level  of  the  photoreceptors  and  retinal  pigment  epithelium 
(RPE)  which  were  not  previously  visible  on  slower,  lower  resolution  systems  (Figure  2.3). 

The  photoreceptor  and  RPE  layers  are  thought  to  be  especially  vulnerable  in  early  age-related 
macular  degeneration,  which  is  a  major  cause  of  vision  loss  in  the  developed  world.  The  features 
in  the  outer  retina  must  be  quantitatively  characterized  in  order  to  provide  establish  a  definition 
of  the  “normal  eye”.  A  nonnative  baseline  is  important  because  it  enables  more  accurate 
assessment  of  subtle  variations  from  the  nonn,  which  may  be  early  indicators  of  disease.  The 
eventual  goal  of  this  work  is  to  identify  quantitative  parameters  that  may  be  used  to  distinguish 
between  the  nonnal  retina  and  early  disease. 
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Figure  2.3.  High  speed  UHR-OCT  imaging  of  the  normal  retina.  New  features  at 
the  level  of  the  photoreceptor  layer  and  retinal  pigment  epithelium  are  visible,  (a) 
Cone  dominated  foveal  region  of  the  eye.  (b)  Rod  dominated  area  of  the  retina. 
ELM  external  limiting  membrane,  IS/OS  inner/outer  photoreceptor  junction,  COS 
cone  outer  segments,  ROS  rod  outer  segments,  RPE  retinal  pigment  epithelium, 
BM  Bruch’s  membrane. 


The  distribution  of  thicknesses  of  different  outer  retinal  layers  in  the  macular  region  in  a  group  of 
48  eyes  of  3 1  normal  human  subjects  in  the  age  range  between  19  and  80  years  was  investigated. 
For  segmentation  of  the  outer  photoreceptor  layer,  retinal  pigment  epithelium  layer  (RPE)  as 
well  as  the  whole  outer  retinal  complex  we  used  scanning  protocol  consisting  of  24  radial  scans 
(Figure  2.3  c)  and  segmentation  software  developed  in  our  group.  This  software  allowed  us  to 
perform  thickness  maps  and  quantitatively  investigate  correlation  of  different  layers  thicknesses 
with  age.  Example  results  of  this  study  are  shown  in  Figure  2.4. 

2.3  Ultrahigh  resolution  OCT  study  of  retinal  pathologies 

UHR-OCT  provides  new  information  about  retinal  structure  that  complements  standard 
diagnostic  techniques  such  as  fluorescein  angiography  (FA),  indocyanine  green  angiography 
(ICG)  or  fundus  photography  by  providing  cross-sectional  images  of  micron-scale  pathological 
changes.  To  take  full  advantage  of  this  ability  in  ophthalmology,  systematic  cross-sectional  and 
longitudinal  studies  of  pathologies  are  necessary.  These  studies  serve  to  illustrate  features 
characteristic  of  various  retinal  pathologies  and  aid  ophthalmologists  in  the  interpretation  of 
OCT  images.  Detailed  qualitative  and  quantitative  analysis  of  UHR-OCT  images  also  enables 
characterization  of  early  markers  of  pathologies  and  improved  understanding  of  pathogenesis. 
Our  group  is  conducting  clinical  investigations  of  different  retinal  disease  in  collaboration  with 
New  England  Eye  Center  (NEEC)  of  Tufts-New  England  Medical  Center  (NEMC)  and 
University  of  Pittsburgh  Medical  Center  (UPMC).  The  study  was  approved  by  the  IRB 
committees  of  MIT,  NEMC  and  UPMC  and  is  compliant  with  the  Health  Insurance  Portability 
and  Accountability  Act  of  1996. 
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Figure  2.4.  Study  on  normal  subjects,  (a)  Selected  macular  regions  within  which 
quantitative  analysis  was  performed:  A  =  0.75  mm,  B=1.5  mm,  C=2.25  mm,  D=3.0 
mm.  (b)  Segmentation  of  outer  retinal  layers:  1  -  inner/outer  photoreceptor 
junction,  2  -  retinal  pigment  epithelium,  3  -  photoreceptor  outer  segment,  4  -  outer 
retinal  complex,  (c)  Example  thickness  map  of  the  photoreceptor  outer  segment  (3) 
of  normal  eye.  (d)  Graph  showing  correlation  within  age  of  retinal  pigment 
epithelium  layer  (2)  thickness  in  the  foveal  region  of  the  macula,  (e)  Table  of 
quantitative  measurements  (mean  values  and  standard  deviations)  of  retinal  outer 
complex  (4)  thickness  in  different  macular  regions. 


A  cross  section  of  retinal  diseases  including:  macular  holes  [72,  75-79],  glaucoma  [80-83],  age- 
related  macular  degeneration  [72,  78,  79,  84],  epiretinal  membranes  [72,  78,  79],  retinal 
dystrophies  [85,  86],  photoreceptor  atrophies  [87],  retinal  detachments  [72,  88],  vitreous 
detachments  [78,  89]  central  serous  chorioretinopathy  [78,  79]  retninis  pigmentosa  [72,  78,  90], 
choroidal  neovascularization  [78,  79],  and  others  have  been  investigated  to  date. 

2.4  Data  visualization  and  analysis 

Ultrahigh  resolution  OCT  data  should  be  visualized  in  a  way  that  reveals  all  information  relevant 
to  diagnosis  of  retinal  pathologies.  In  the  simplest  visualization  method  single  B-scans  (or  OCT 
frames)  showing  the  most  representative  cross-sections  of  investigated  pathology  are  presented 
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using  an  appropriate  color  scale  (commonly  gray  or  false  color  scales  representing  different 
levels  of  light  scattering  in  eye  structures).  Such  a  solution  is  usually  utilized  in  standard 
resolution  StratusOCT  or  ultrahigh  resolution  time  domain  OCT,  where  image  acquisition  speeds 
are  relatively  slow.  However,  with  the  advent  of  high-speed  imaging,  where  as  many  as  ~200 
images  are  acquired  at  once,  more  advanced  methods  of  visualization  and  data  display  are 
required. 

Improvements  in  imaging  speed  enabled  by  spectral/Fourier  domain  OCT  allow  raster  scan 
protocols  (Figure  2.2  b)  to  obtain  densely  sampled  three-dimensional  OCT  data.  This  data  can  be 
visualized  using  a  range  of  novel  visualization  methods  [71,  79].  For  example,  it  can  be  used  to 
create  OCT  fundus  images  by  summing  the  three-dimensional  data  set  along  the  axial  direction 
at  each  transverse  {en  face)  position  on  the  retina  (Figure  2.5.  A,  panel  a).  This  corresponds  to 
detecting  ah  of  the  light  backscattered  or  backreflected  from  each  fundus  position,  and  therefore 
is  analogous  to  fundus  photography  or  scanning  laser  ophthalmoscopy.  Since  the  OCT  fundus 
image  is  generated  directly  from  the  3-D  OCT  data,  OCT  images  are  precisely  and  reproducibly 
registered  with  the  fundus.  The  result  is  then  displayed  as  a  grey  scale  image  that  enables  direct 
comparison  of  OCT  findings  with  those  from  clinical  examination,  such  as  fundus  photographs 
or  fluorescein  angiography.  It  is  also  possible  to  generate  OCT  fundus  images  selectively 
displaying  specific  retinal  layers  or  specific  retinal  features  (Figure  2.5.  A,  panel  b). 

Another  set  of  different  visualization  methods  is  shown  in  Figure  2.5.B.  Before  visualization, 
individual  frames  from  the  three-dimensional  data  set  were  correlated  to  remove  axial  motion 
artifacts.  Correlation  between  consecutive  axial  scans  is  unnecessary  due  to  the  fast  acquisition 
rate.  Virtual  cross-sectional  images,  segmented,  virtual  perspective,  and  cut-away  rendered  views 
of  the  intraretinal  layers  were  in  this  case  generated  using  volumetric  rendering  software  similar 
to  that  for  magnetic  resonance  image  processing. 

Three  dimensional  data  sets  of  radial  scans  densely  covering  the  retina  (Figure  2.2  c)  may  be 
used  for  quantitative  analysis  of  OCT  tomograms  [71,  82,  91,  92],  Segmentation  of  different 
structures  allows  for  calculation  of  volume  of  selected  retinal  features  (Figure  2.5.B,  panel  f)  as 
well  as  mapping  intraretinal  layer  thicknesses.  The  later  is  especially  important  in  quantifying 
macular  edema  which  is  a  consequence  of  many  pathological  retinal  conditions  such  as  diabetic 
retinopathy,  epiretinal  membrane  formation,  ocular  inflammation,  retinal  vascular  occlusion,  and 
cataract  extraction.  Macular  thickness  analysis  is  particularly  helpful  in  guiding  and  assessing 
treatments  for  macular  edema  such  as  photodynamic  therapy,  intravitreal  corticosteroids,  and 
vitrectomy.  It  is  also  important  for  the  detection  and  monitoring  of  glaucoma  [82]. 

An  example  of  retinal  thickness  analysis  is  shown  in  Figure  2.6.  The  patient  had  resolving 
central  serous  chorioretinopathy  (CSCR)  in  the  left  eye.  No  pathological  changes  were  observed 
in  the  right  eye.  The  patient  was  imaged  with  high-speed  UHR  OCT  using  24  radial  scans 
(Figure  2.2  c).  Figures  2.6  a  show  horizontal  6  mm  OCT  image  through  the  macular  region.  The 
thickness  of  ah  retinal  layers,  outer  retinal  complex  as  well  as  retinal  pigment  epithelium  (RPE) 
was  mapped.  The  thickness  maps  of  the  RPE  of  both  eyes  are  shown  in  Figures  2.6  c.  Thinning 
of  the  RPE  in  the  left  eye  is  clearly  visible. 
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Figure  2.5.  Visualization  in  high  speed,  UHR  OCT  imaging.  (A)  Full  thickness 
macular  hole:  a.  OCT  fundus  image  created  by  summation  data  in  the  axial 
direction;  b.  Coronal  section  through  macular  hole  (red  line  in  Figure  c), 
demonstrating  “mosaic-like”  cystic  pattern  surrounding  macular  hole;  c.  Cross- 
sectional  view  of  macular  hole.  B.  Volumetric  analysis  of  high-speed,  UFIR  OCT 
data  of  full  thickness  macular  hole  with  epiretinal  membrane,  a.  Volume 
rendering  of  the  retina;  b.  Visualization  of  segmented  epiretinal  membrane  along 
with  the  rendered  volume;  c.  OCT  image  chosen  from  the  three  dimensional 
dataset  showing  epiretinal  membrane;  d.  Segmentation  of  the  epiretinal  membrane 
in  the  same  image;  e.  Map  of  epiretinal  membrane  registered  to  the  OCT  fundus 
image;  f.  Segmentation  of  macular  hole  and  calculation  of  hole  volume. 
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Figure  2.6.  Macular  thickness  analysis.  Example  of  a  patient  with  resolving 
central  serous  chorioretinopathy  in  the  left  eye  (B),  and  healthy  right  eye  (A),  (a) 
One  of  the  OCT  cross  sections  (horizontal)  taken  from  the  set  of  24  radial  scans 
(refer  to  the  Figure  2.2  c).  (b)  Segmentation  of  different  retinal  layers  in  the 
macular  region,  (c)  Thickness  maps  of  the  retinal  pigment  epithelium  (RPE): 
thinning  of  the  RPE  in  the  left  eye  is  visualized. 


2.5  Small  animal  retinal  imaging 

Research  on  ocular  diseases  is  limited  by  the  restrictions  on  studying  pathophysiologic  processes 
in  the  human  eye.  Murine  (rat  and  mouse)  models  of  ocular  disease  therefore  provide  powerful 
tools  for  characterization  of  disease  pathogenesis  and  response  to  treatment.  While  enucleation 
and  histology  are  the  gold  standard  for  characterization  of  microstructural  changes  in  animals, 
non-invasive  structural  imaging  has  the  potential  to  reduce  the  need  for  sacrifice  and  histology  in 
many  studies.  New  OCT  detection  methods  enable  high-speed  imaging  as  well  as  numerical 
dispersion  compensation  and  spectral  shaping,  which  are  crucial  for  visualizing  the  fine 
structures  in  the  retina. 

A  high-speed,  UHR-OCT  system  was  developed  for  small  animal  retinal  imaging.  This  system 
achieves  imaging  speeds  of  24,000  axial  scans  per  second,  an  improvement  of  ~100x  over 
previous  UHR-OCT  systems.  Using  a  broadband  superluminescent  diode  light  source,  an  axial 
image  resolution  of  2.8  pm  is  achieved.  A  10  pm  transverse  resolution  in  air  is  obtained  using  a 
post-objective  scanning  microscope.  High-speed  UHR-OCT  enables  high  quality  imaging  of  the 
murine  retina  and  the  visualization  of  all  major  intraretinal  layers.  Raster  scan  protocols  enable 
volumetric  imaging  and  comprehensive  coverage  of  a  region  of  the  retina.  An  OCT  fundus 
image,  akin  to  a  fundus  photograph,  can  be  generated  by  axial  summation  of  three-dimensional 
OCT  data,  and  enables  precise  registration  of  OCT  measurements  to  retinal  fundus  features. 
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Figure  2.7  shows  an  example  of  an  OCT  cross  sectional  image  (A)  and  a  fundus  image  (B), 
which  shows  features  such  as  retinal  blood  vessels.  Three-dimensional  imaging  enables 
quantification  of  retinal  structure,  which  promises  to  allow  repeated,  non-invasive  measurements 
to  track  disease  progression,  reducing  the  need  for  sacrifice  and  histology.  This  capability  can 
accelerate  the  translation  from  basic  research  studies  in  rats  and  mice  into  clinical  care. 


Figure  2.7.  An  OCT  cross-sectional  image  of  the  rat  retina  (A)  with  the  OCT 
fundus  image  (B),  obtained  by  summation  of  3D-OCT  data.  The  location  of  the 
cross-sectional  image  on  the  OCT  fundus  image  is  shown  by  a  white  line. 

2.6  Functional  Retinal  Imaging 

Functional  impairment  may  precede  structural  changes  in  a  number  of  retinal  diseases  including 
age-related  macular  degeneration,  glaucoma,  and  diabetic  retinopathy.  Objective  methods  for 
detecting  functional  impairment  in  the  retina  can  improve  sensitivity  and  reproducibility  for 
tracking  of  changes  over  time.  There  are  several  objective  methods  for  measuring  the  functional 
response  of  the  retina  and  optic  nerve  including  multifocal  electroretinography,  pattern 
electroretiography,  multifocal  visual  evoked  potentials,  and  pupil  perimetry. 

Recently,  optical  techniques  have  been  demonstrated  to  provide  physiological  infonnation.  Small 
changes  in  optical  properties  caused  by  membrane  depolarization  and  cell  swelling  have  been 
measured  in  the  cortex  of  the  brain  [93,  94].  Changes  in  oxygen  saturation  have  been  measured 
in  response  to  retinal  stimulation  using  fundus  reflectometry  [95].  Other  studies  demonstrated 
changes  in  retinal  reflectance  in  response  to  a  flash  stimulus  that  correlated  with  multifocal  ERG 
[96].  While  these  fundus  imaging  techniques  provide  good  spatial  resolution  in  the  transverse 
direction,  they  do  not  have  depth  resolving  capability. 

OCT  imaging  of  neural  activity  was  demonstrated  in  the  sea  slug  ganglion  [97]  as  well  as  in  the 
cat  visual  cortex  [98,  99].  Recent  work  using  conventional  OCT  demonstrated  measurement  of 
retinal  functional  changes  in  response  to  a  flash  stimulus  [100,  101].  These  studies  showed 
changes  in  OCT  signals  in  the  photoreceptor  inner  and  outer  segments  and  were  performed  using 
ex  vivo  preparations.  The  investigators  attributed  the  observed  changes  to  photoreceptor  disc 
membrane  hyperpolarization  in  response  to  a  strong  stimulus,  although  this  hypothesis  remains 
to  be  tested. 
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Using  high-speed,  ultrahigh  resolution  OCT  with  2.8  pm  resolution  and  24,000  axial  scans  per 
second,  we  demonstrated  in  vivo  measurement  of  functional  reflectance  changes  in  the  rat  retina. 
The  rat  retinal  model  was  chosen  because  of  its  widespread  use  in  basic  research  on  ocular 
disease.  We  have  developed  a  novel  measurement  protocol  that  enables  the  simultaneous 
correction  of  motion  artifacts  and  speckle  averaging.  This  protocol  enables  a  dramatic  reduction 
of  physiologic  noise  from  motion  and  allows  single  shot  measurements  of  photoreceptor 
function,  without  the  need  to  average  across  multiple  trials.  The  spatial  and  spectral  dependence 
of  the  functional  reflectance  change  was  characterized,  as  well  as  the  effect  of  photopigment 
bleaching. 

For  functional  OCT  experiments,  animals  are  dark-adapted  for  12  hours,  anesthetized,  and 
prepared  under  dim  red  light.  A  transient  white  light  stimulus  is  delivered  during  the  functional 
recording,  while  a  selected  region  of  the  retina  is  repeatedly  monitored  for  reflectance  changes. 
The  spectrum  of  wavelengths  used  in  the  OCT  measurements  was  from  817-962  mn,  well 
outside  the  sensitivity  range  of  the  rat  retina.  The  stimulus  duration  and  intensity  were  chosen  to 
prevent  full  photoreceptor  bleaching.  Figure  2.8  shows  the  amplitude  reflectance  from  the 
photoreceptor  outer  segments,  plotted  as  a  function  of  time.  As  shown,  the  stimulus  induces  a 
— 10-15%  increase  in  the  average  amplitude  reflectance  from  the  photoreceptor  outer  segments. 

OCT  may  eventually  become  a  clinical  tool  for  assessing  photoreceptor  health  in  the  early  stages 
of  retinal  disease.  We  have  shown  that  single  shot  measurements  can  be  achieved  in  the 
anesthetized  rat.  Our  results  and  methodology  are  promising  for  the  eventual  application  of  these 
techniques  in  the  ophthalmology  clinic.  OCT  structural  imaging  can  be  performed  to 
longitudinally  follow  areas  of  functional  deficit  for  structural  changes,  or  vice  versa.  The 
characterization  of  the  response  topography  and  temporal  dependence  in  the  rat  model  will  be 
useful  in  developing  protocols  for  functional  OCT  measurements  in  humans.  Although  the 
immediate  application  of  functional  OCT  measurement  technologies  to  humans  is  attractive, 
further  systematic  in  vivo  studies  in  a  well-controlled  retinal  model  such  as  the  rat  retina  are 
important  to  characterize  the  observed  optical  changes. 


atimotua  duration 

095 1 _ x _ i _ wp— MfM-s _ i _ i _ . _ , _ x_ 

0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 


time  (s) 

Figure  2.8.  A  light  stimulus  induces  an  increase  in  the  average  amplitude 
reflectance  from  the  photoreceptor  outer  segments.  Each  data  point  in  this  figure 
corresponds  to  the  average  reflectance  from  the  photoreceptor  outer  segments 
obtained  from  a  small  region  of  the  retina. 
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3.  Endoscopic  OCT  Studies 


3.1  Endoscopic  Imaging  of  Barrett’s  Esophagus 

Optical  coherence  tomography  (OCT)  can  generate  high  resolution,  cross-sectional  images  of 
biological  tissues  in  situ  and  in  real  time  [1,2,  102].  OCT  can  function  as  a  type  of  optical  biopsy 
to  enable  imaging  of  tissue  microstructure  with  the  resolution  approaching  that  of  standard 
excision  biopsy,  without  the  need  of  excising  the  tissue  specimen  [3-5].  One  promising 
application  of  optical  biopsy  using  OCT  is  the  endoscopic  imaging  of  the  gastrointestinal  (GI) 
tract.  In  contrast  to  conventional  endoscopy,  which  can  only  visualize  the  surface  alterations, 
OCT  can  detect  changes  in  tissue  morphology  beneath  the  tissue  surface.  Therefore,  endoscopic 
imaging  with  high  resolution  OCT  could  potentially  improve  the  detection,  visualization,  and 
diagnosis  of  gastrointestinal  diseases. 

Several  researchers  have  investigated  endoscopic  OCT  imaging  in  the  human  gastrointestinal 
(GI)  tract,  including  the  esophagus  and  stomach,  the  small  and  large  intestine,  and  the  bile  duct 
[103-113],  In  addition  to  structural  imaging,  imaging  of  blood  flow  can  also  be  performed  using 
endoscopic  Doppler  OCT  [114].  OCT  imaging  has  demonstrated  high  sensitivities  and 
specificities  (100%  sensitivity  and  93%  specificity  in  a  retrospective  study,  97%  sensitivity  and 
92%  specificity  in  a  prospective  study)  for  the  identification  of  specialized  intestinal  metaplasia 

[115] ,  OCT  imaging  could  also  distinguish  hyperplastic  from  adenomatous  polyps  in  the  colon 

[116] ,  Recent  studies  investigated  the  capability  of  OCT  for  detecting  dysplasia  in  Barrett’s 
esophagus.  Evans  et  al.  reported  a  sensitivity  of  83%  and  specificity  of  75%  for  detection  of 
dysplasia  with  blinded  scoring  of  OCT  images  [117].  Isenberg  et  al.  reported  an  accuracy  of  78% 
for  the  detection  of  dysplasia  in  patients  with  Barrett's  esophagus  [118].  However,  currently 
almost  all  clinical  studies  have  been  perfonned  using  standard  OCT  with  10-15  pm  resolution. 
Ultrahigh  resolution  OCT  could  enhance  the  imaging  perfonnance  for  the  identification  of  early 
neoplastic  changes,  and  could  improve  the  sensitivity  of  biopsy  by  reducing  false  negative  rates 
from  sampling  errors. 

Our  group  developed  a  research  prototype  ultraghigh  resolution  endoscopic  OCT  system  for  use 
in  the  endoscopy  suite.  In  order  to  achieve  the  high  powers  and  short  coherence  lengths 
necessary  for  high  resolution,  high  speed  imaging,  a  Cr:Forsterite  laser  was  used  as  the  light 
source  [25].  This  laser  generates  broadband  spectrum  in  the  1300  nm  wavelength  regime.  The 
output  bandwidth  can  be  increased  by  using  nonlinear  effects  in  optical  fibers  to  yield  a 
coherence  length  of  5  pm  or  less.  Due  to  the  bandwidth  limitations  in  the  optical  components  in 
the  sample  and  reference  arms,  the  back-coupled  spectrum  on  the  detector  had  a  bandwidth  of 
150  nm,  which  corresponds  to  a  theoretical  axial  resolution  of  4.6  pm  in  air.  The  width  of  the 
measured  axial  point  spread  function  is  5  pm,  corresponding  to  ~4  pm  resolution  in  the  tissue, 
which  is  2-3x  finer  than  standard  OCT  systems.  The  system  sensitivity  was  102  dB  at  4  Hz 
frame  rate  with  up  to  15  mW  power  on  the  sample. 

Working  in  collaboration  with  Dr.  Hiroshi  Mashimo,  MD,  PhD  at  the  Harvard  Medical  School, 
we  performed  a  study  to  investigate  endoscopic  OCT  imaging.  Ultrahigh-resolution  (UHR)  OCT 
imaging  was  performed  in  the  esophagus  and  stomach  of  patients  undergoing  upper  GI 
endoscopy  at  the  Boston  Veteran  Affairs  Medical  Center  (VAMC).  Informed  consent  was 
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obtained  from  patients  enrolled  under  a  protocol  approved  jointly  by  the  Institutional  Review 
Board  (IRB)  of  Boston  VAMC,  the  Committee  on  Human  Studies  at  Harvard  Medical  School 
and  the  Committee  on  the  Use  of  Humans  as  Experimental  Subjects  (COUHES)  of 
Massachusetts  Institute  of  Technology.  Imaging  was  performed  using  a  1.8  mm  OCT  imaging 
catheter  which  was  introduced  through  the  accessory  channel  of  the  endoscope.  The  OCT 
imaging  catheter  was  disinfected  prior  to  use.  While  in  the  field  of  view  of  the  endoscope,  the 
OCT  imaging  catheter  was  positioned  to  the  areas  of  interest  by  maneuvering  the  tip  of  the 
endoscope. 

A  total  of  50  patients  with  previous  history  of  Barrett’s  esophagus  were  imaged.  Images  from 
both  normal  and  abnormal  areas  in  the  upper  GI  tract  were  acquired  and  later  classified 
according  to  the  results  of  pinch  biopsy  histology.  The  pathologic  diagnoses  from  biopsy 
specimens  included:  Barrett’s  esophagus  without  dysplasia  (17),  indefinite  for  dysplasia  (5), 
low-grade  dysplasia  (16),  high-grade  dysplasia  (4)  and  adenocarcinoma  (2). 


Figure  3.1.  A,  B:  In  vivo  images  of  normal  esophagus  acquired  using  ultrahigh 
resolution  endoscopic  OCT  (wavelength  1.3  pm,  resolution  5  pm  axial  x  15  pm 
transverse).  Layered  architecture  of  the  esophageal  wall  including  the  epithelium 
(ep),  lamina  propria  (lp),  muscularis  mucosa  (mm),  submucosa  (sm),  muscularis 
propria  (mp)  are  visualized. 

Figure  3.1  shows  UHR  OCT  images  of  normal  esophagus  with  squamous  epithelium.  Figure 
3.1  A  reveals  the  characteristic  five-layer  architecture  of  normal  esophagus  mucosa  which  has 
been  observed  in  previous  standard  resolution  OCT  imaging  studies.  The  UHR  OCT  images 
show  the  relatively  homogeneous  epithelium  (ep),  the  high-backscattering  band  (brighter)  of  the 
lamina  propria  (lp),  the  low-backscattering  (darker)  muscularis  mucosae  (mm),  the  high- 
backscattering  submucosa  (sm),  and  the  low-backscattering  and  thick  muscularis  propria  (mp).  A 
three-layered  architecture  was  also  identified  in  the  same  patient  as  shown  in  Figure  3. IB.  The 
high-backscattering  connective  tissue  layers  of  the  lamina  propria  and  submucosa  merge  without 
the  lower  back-scattering  muscularis  mucosae  between  them.  This  finding  is  consistent  with  the 
fact  that  the  longitudinal  fibers  of  the  muscularis  layer  occur  in  bundles,  and  a  given  cross- 
sectional  OCT  image  may  not  intercept  one  of  these  bundles.  The  identification  of  the  muscularis 
mucosa  is  important  for  staging  of  esophageal  cancer.  Although  these  features  have  been 
observed  in  standard  resolution  OCT,  they  are  sharper  and  easier  to  identify  in  UHR  OCT 
images. 

Representative  endoscopic  images  with  the  catheter  placement  and  resulting  UHR  OCT  images 
of  regions  of  Barrett’s  esophagus  from  different  patients  are  presented  in  Figure  3.2.  The  UHR 
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OCT  images  demonstrated  clear  differences  in  the  tissue  architecture  of  the  Barrett's  region  when 
compared  to  normal  esophageal  squamous  mucosa.  The  horizontally  layered  esophageal 
squamous  epithelium  was  replaced  by  the  presence  of  crypt-like  glandular  architectures,  which 
are  more  heterogeneously  scattering  compared  to  normal  esophagus.  Low-backscattering  glands 
were  frequently  observed  within  Barrett’s  mucosa,  with  interlaced  regions  of  high-backscattering 
connective  tissue  corresponding  to  the  lamina  propria.  Although  individual  patients  exhibited 
variation  in  the  detailed  architecture  of  Barrett’s  metaplasia  as  visualized  with  UHR  OCT,  the 
common  OCT  image  features  of  crypt-like  glandular  architectures  were  present  in  all  patients 
with  Barrett’s  esophagus.  The  deeper  horizontal  stratified  muscularis  mucosae  and  submucosal 
layers  of  the  normal  esophagus  were  preserved  in  patients  with  Barrett’s  esophagus. 


Figure  3.2.  Ultrahigh  resolution  endoscopic  OCT  imaging  of  Barrett’s  esophagus 
from  2  patients  showing  variation  in  Barrett’s  morphology:  A,  C:  Endoscopic 
view;  B,  D:  OCT  images  (wavelength  1.3  pm,  resolution  5  pm  axial  x  15  pm 
transverse).  A,  B  are  from  the  same  subject;  C,  D  are  from  the  same  subject. 

Figure  3.3  shows  a  case  of  esophageal  adenocarcinoma  underneath  benign  squamous  epithelium 
in  a  patient  status  post  photodynamic  therapy  (PDT).  The  irregularly  backscattering  neoplastic 
mucosa  of  the  adenocacinoma  is  dramatically  different  from  that  of  normal  and  Barrett’s 
esophagus.  OCT  images  of  adenocarcinoma  show  heterogeneous  scattering  corresponding  to 
prominent  morphologic  disorganization,  irregularly  shaped  and  crowded  glandular  architecture, 
with  reduced  overall  back-scattering  intensity.  Features  of  esophageal  adenocarcimona 
visualized  with  UHR  OCT  agree  with  those  found  in  previously  published  OCT  imaging  studies 
and  are  consistent  with  pathological  findings.  Qualitatively,  the  OCT  images  show  progressive 
increase  of  architectural  irregularity  from  Barrett’s  esophagus  to  high-grade  dysplasia  and 
eventually  to  adenocarcinoma.  A  comparison  was  also  performed  to  investigate  imaging  results 
between  endoscopic  ultrasound  (EUS)  and  UHR  OCT.  Figure  3.3  shows  both  types  of  images 
derived  from  the  same  patient.  UHR  OCT  demonstrates  higher  imaging  resolution  for  epithelial 
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architecture  than  standard  radial  EUS.  However,  EUS  has  higher  image  penetration  than  OCT, 
enabling  the  depth  of  the  tumor  invasion  to  be  assessed.  This  result  suggests  that  these  two 
imaging  modalities  might  provide  complementary  information  for  the  detection  and  staging  of 
esophageal  adenocarcinoma. 


Figure  3.3.  Ultrahigh  resolution  endoscopic  OCT  imaging  of  adenocarcinoma 
from  a  patient  status  post  PDT.  A,  Endoscopic  view;  B,  OCT  image  (wavelength 
1.3  pm,  resolution  5  pm  axial  x  15  pm  transverse);  C,  Endoscopic  ultrasound 
images  of  the  same  area;  D,  Corresponding  histology  (H&E,  ori.  mag.  x40). 

These  studies  demonstrate  the  ability  of  UHR  OCT  to  differentiate  normal  and  pathologic  tissue. 
The  improved  image  resolution  of  OCT  enabled  the  visualization  of  architectural  morphology 
features  such  as  the  normal  layered  structure  of  the  epithelium  versus  glandular  and  columnar 
structures  associated  with  Barrett’s  esophagus  more  clearly.  Further  and  more  detailed 
investigations  using  ultrahigh  resolution  OCT  imaging  in  clinical  endoscopic  studies  will  be 
necessary  to  evaluate  the  ultimate  impact  of  improved  resolution  on  the  diagnostic  capabilities 
for  detecting  Barrett’s  esophagus  associated  high-grade  dysplasia  and  carcinoma.  In  the  context 
of  surveillance  of  patients  with  Barrett’s  for  high  grade  dysplasia  and  adenocarcinoma,  the  most 
intriguing  application  of  OCT  would  be  to  direct  excisional  biopsy  to  reduce  sampling  errors. 
One  can  envision  new  OCT  imaging  probes  which  integrate  OCT  imaging  with  pinch  biopsy  to 
provide  a  real  time  “first  look”  at  pathology  prior  to  excision  and  processing  of  a  specimen.  If 
successful,  this  technology  could  be  used  to  help  guide  excision  biopsy  to  reduce  sampling  errors 
and  false  negative  rates.  This  could  improve  sensitivity  of  diagnosis,  reduce  the  cost  of 
surveillance  and  provide  enhanced  diagnosis  and  treatment  decisions. 
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3.2.  Image  Analysis  Methods 


Image  analysis  has  been  extensively  applied  in  various  biomedical  imaging  technologies,  such  as 
ultrasound,  magnetic  resonance  imaging  (MRI),  and  computerized  tomography  (CT),  to  assist  in 
the  diagnosis  of  diseases.  Previous  studies  have  demonstrated  the  diagnostic  capability  of 
endoscopic  OCT  in  Barrett’s  esophagus  [115],  esophageal  dysplasia  [117,  118]  and  colon 
dysplasia  [116]  by  assessment  of  tissue  morphological  disorganization.  However,  such  human- 
based  diagnosis  approaches  would  require  the  special  training  in  interpreting  the  OCT  images, 
might  result  in  the  intra-  and  inter-observer  variations,  and  potentially  become  time-consuming  if 
the  number  of  images  is  large.  Recently,  computer-based  image  analysis  approaches  have  been 
applied  to  OCT  images  for  classification  and  diagnosis  [119,  120].  Preliminary  studies  have 
demonstrated  the  promise  of  diagnosis  of  dysplasia  from  Barrett’s  esophagus  using  texture 
analysis  on  standard  resolution  OCT  images  (15  pm  axial  resolution)  [120],  This  approach  will 
potentially  minimize  the  intra-  and  inter-observer  variations,  and  enable  automatic  analysis  of 
large  volume  of  data  sets  through  computer  aided  diagnosis  (CAD). 

We  have  investigated  the  capability  of  automatic  tissue  characterization  and  diagnosis  using 
texture  analysis  for  both  standard-resolution  (SDR-OCT)  and  ultrahigh-resolution  OCT  (UHR- 
OCT)  in  detection  of  Barrett’s  esophagus  (BE).  Image  features  such  as  coarseness,  contrast, 
periodicity,  and  orientation  can  be  quantified  using  texture  analysis.  Features  of  50  endoscopic 
OCT  images  (25  nonnal  and  25  Barrett’s  esophagus)  were  extracted  through  computational 
metrics  including  the  spatial  frequency  analysis  using  two-dimensional  Fourier  transform  (2D 
FFT),  and  statistical  texture  analysis  using  center  symmetric  auto-correlation  [121]  (see  Equation 
3.1).  These  features  were  further  processed  using  principle  component  analysis  (PCA)  to  reduce 
the  variable  dimensions  and  increase  the  discriminative  power. 


SCOV  =  1/41  (gru)(gl’-u) 

VAR  =1/81  (gf-gi’2)-u2 
BVAR  =  1/161  (gi+gi’)2-u2 
WVAR=  1/161  (grgi’)2 
SVR  =  WVAR/BVAR 
SAC  =  SCOV/VAR 
*u:  local  mean  value 

Equation  3.1.  Algorithm  for  calculating  the  image  features  through  center 
symmetric  auto-correlation.  SCOV  gray  scale  texture  covariance,  SAC 
nonnalized  SCOV,  BVAR  between-pair  variance,  WVAR  within-pair  variance, 
VAR:  local  variance,  SVR  variance  ratio. 
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Figure  3.4  plots  the  first  two  principle  components  (T 1  and  T2)  from  PCA  for  normal  and 
Barrett’s  esophagus  image  features.  T1  and  T2  are  linear  combinations  of  original  variables, 
including  2D  FFT  features  and  statistical  texture  features.  The  quantitative  metrics  are:  1)  the 
variable  reduction;  2)  separation  capability.  Variable  reduction  is  estimated  by  the  variances  of 
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first  two  principle  components  (T 1  and  T2).  UHR-OCT  shows  more  distinct  difference  in  T1  and 
T2,  suggesting  data  are  more  spread  out  in  one  dimension  therefore  facilitates  the  further 
discrimination  (see  Table  3.1).  The  discrimination  capability  for  normal  and  BE  groups  is 
estimated  by  t-statistic  which  is  the  ratio  of  the  distance  between  the  centroid  of  two  groups 
(between-group)  divided  by  the  weighted  standard  error  of  each  group  (within- group).  UHR 
shows  larger  ratio  of  between-group  variance  to  within-group  variance  than  SDR  (Table  3.2), 
suggesting  the  normal  and  Barrett’s  esophagus  data  are  further  separated  in  UHR  OCT  images. 
These  results  suggest  that  the  improved  axial  resolution  will  enhance  the  capability  in  discerning 
the  fine  features  in  OCT  images,  which  could  be  useful  in  diseases  diagnosis.  The  enhanced 
classification  abilities  using  UHR  OCT  would  help  in  the  computer  aided  diagnosis  of  GI 
diseases.  This  work  demonstrates  enhanced  discrimination  of  Barrett’s  esophagus  from  normal 
esophagus,  but  these  results  may  also  prove  important  in  the  discrimination  of  dysplasia  from 
Barrett’s  tissues.  Similar  analysis  algorithms  have  been  applied  to  high-grade  dysplasia  images 
and  the  analysis  is  undergoing. 
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Figure  3.4.  Scatter  plots  of  first  two  principle  components  (T 1  and  T2)  from  PCA 
for  normal  and  Barrett’s  esophagus.  T1  and  T2  are  linear  combinations  of  original 
variables,  including  2D  FFT  features  (A,B)  and  statistical  texture  features  (C,D). 
A:  Plots  of  T1  and  T2  from  2D  FFT  features  of  SDR  images;  B:  Plots  of  T1  and 
T2  from  2D  FFT  features  of  UHR  images.  The  UHR  features  shows  more 
condensed  along  T1  axis  (more  variance  in  T1  and  less  variance  in  T2  compared 
to  SDR  features,  see  values  in  Table  3.1).  BE  and  normal  groups  are  further 
separated  in  UHR  (see  values  in  Table  3.2).  C:  Plots  of  T1  and  T2  from  texture 
features  of  SDR  images;  D:  Plots  of  T1  and  T2  from  texture  features  of  UHR 
images.  Similarly,  the  UHR  features  shows  more  concentrated  along  T1  axis  and 
further  separated  for  BE  and  normal  groups. 
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Table  3.1.  Comparison  of  Variable  Reduction 


Variance 

SDR 

(FFT) 

UHR  (FFT) 

SDR 

(Texture) 

UHR 

(Texture) 

T1 

3.02 

4.23 

5.57 

7.21 

T2 

2.23 

0.88 

2.10 

1.06 

Table  3.2.  Comparison  of  Discrimination 


t- Statistics 

SDR 

UHR 

FFT 

7.07 

16.05 

Texture 

10.40 

37.77 

4.  Optical  Coherence  Microscopy 
4.1  Introduction 

Ultrahigh  resolution  Optical  Coherence  Tomography  (OCT)  can  achieve  1-2  pm  axial  resolution 
in  tissue,  but  is  limited  in  transverse  resolution  due  to  the  low  numerical  aperture  (NA)  focusing 
to  maintain  a  sufficient  depth  of  field  over  the  range  of  the  cross-sectional  image  [8],  The 
relatively  low  lateral  resolution  achievable  with  cross-sectional  OCT  is  generally  insufficient  for 
imaging  of  cellular  features  and  therefore  limits  the  utility  of  OCT  in  applications  requiring 
cellular  level  diagnostics.  To  extend  the  imaging  power  of  OCT  to  very  high  transverse 
resolution,  we  are  developing  a  technology  known  as  Optical  Coherence  Microscopy  (OCM), 
which  combines  OCT  with  confocal  microscopy.  Figure  4.1  compares  the  image  penetration 
depth  and  resolution  of  OCM  with  that  of  OCT  and  confocal  microscopy.  OCM  can  provide 
enhanced  penetration  depth  compared  to  standard  confocal  microscopy  while  dramatically 
improving  the  resolution  over  typical  cross-sectional  OCT  imaging  methods. 

Optical  Coherence  Microscopy  overcomes  the  depth  of  field  limitation  present  in  traditional 
OCT  imaging  by  imaging  in  the  en  face  plane  rather  than  the  cross-sectional  plane.  To  image  en 
face,  the  optical  path  length  of  the  reference  arm  is  matched  exactly  to  the  focus  of  the  sample 
arm  microscope  while  scanning  a  transverse  raster  pattern  on  the  tissue.  This  eliminates  the  need 
for  path  length  scanning  to  generate  an  axial  depth  map  and  allows  the  use  of  high  NA  lenses  to 
provide  very  small  spot  sizes.  Because  of  the  high  NA  focusing  used  in  OCM,  the  field  of  view 
in  the  images  is  necessarily  smaller  than  in  OCT.  OCT  images  architectural  features  over  a  field 
of  view  of  3-6  mm  square  while  OCM  can  image  down  to  the  cellular  level  but  with  reduced 
field  of  view  of  100-500  pm  square. 

OCM  has  the  unique  advantage  of  using  two  distinct  optical  sectioning  techniques  -  confocal 
gating  and  coherence  gating.  The  multiplicative  effect  of  the  two  sectioning  methods  strengthens 
the  overall  optical  sectioning  power,  allowing  increased  rejection  of  unwanted,  out  of  focus 
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Figure  4.1.  Transverse  resolution  and  image  penetration  in  optical  coherence 
microscopy  (OCM).  OCM  can  enhance  image  penetration  compared  to  confocal 
microscopy  alone  while  significantly  improving  transverse  resolution  in  OCT  to 
enable  cellular  level  imaging. 

scattered  light.  Studies  from  our  group  as  well  as  others  have  demonstrated  that  combined 
confocal  and  coherence  gating  can  provide  improved  imaging  depth  compared  to  confocal  alone 
[122,  123].  The  addition  of  high  sensitivity  coherence  gated  detection  to  confocal  detection 
extends  the  imaging  depth  in  scattering  media  to  the  shot  noise  quantum  limit,  providing  a  factor 
of  2-3  increase  over  standard  confocal  microscopy. 

The  use  of  multiple  optical  sectioning  techniques  also  allows  considerable  flexibility  in  system 
design  for  achieving  high-resolution  cellular  images.  Broad  bandwidth  light  sources  as  used  in 
ultrahigh  resolution  OCT  can  provide  thin  optical  sectioning  via  coherence  gating,  and  the 
confocal  sectioning  can  be  relaxed  to  facilitate  development  of  miniaturized  imaging  devices. 

Figure  4.2  compares  the  confocal  axial  and  transverse  imaging  resolution  as  a  function  of  the 
numerical  aperture  of  the  probe  optics  to  demonstrate  this  operating  limit  for  OCM.  The  axial 
section  thickness  degrades  much  more  quickly  than  the  transverse  resolution,  and  there  exists  a 
region  where  the  transverse  resolution  is  sufficient  for  cellular  imaging  but  the  axial  resolution  is 
not.  Addition  of  a  short  coherence  gate  to  provide  tissue  sectioning  can  therefore  make  cellular 
imaging  possible  with  much  lower  NA  than  is  sufficient  for  confocal  microscopy  alone.  This 
operating  regime  for  OCM  imaging  has  very  important  clinical  implications,  since  it  promises  to 
allow  cellular  imaging  with  small  diameter  probes  compatible  with  standard  endoscopic  and 
laparoscopic  procedures. 

We  have  developed  a  high  speed  OCM  system  capable  of  imaging  with  broadband  femtosecond 
laser  sources  as  used  in  ultrahigh  resolution  OCT.  The  system  used  a  femtosecond  TiiSapphire 
laser  and  a  novel  broadband  optical  phase  modulator  to  provide  high  resolution  imaging  at  800 
mn  wavelength.  Figure  4.3  presents  key  results.  A  short  coherence  gate  of  ~3  pm  was  combined 
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Figure  4.2.  Numerical  aperture  requirement  for  OCM  compared  to  confocal 
microscopy.  OCM  can  image  with  high  transverse  resolution  at  much  lower 
numerical  aperture  than  confocal  microscopy  because  it  does  not  depend  on  high 
axial  resolution  for  optical  sectioning. 


Figure  4.3.  OCM  can  perform  cellular  imaging  with  lower  NA  than  confocal 
microscopy.  (Left)  A  short  coherence  gate  is  used  to  compensate  a  longer  confocal 
gate,  which  will  relax  probe  design  constraints  for  endoscopic  imaging.  (Right)  High 
transverse  resolution  cellular  imaging  can  be  achieved  despite  weaker  confocal 
sectioning.  In  vivo  OCM  cellular  imaging  was  performed  in  Xenopus  laevis  tadpole 
(A)  and  human  skin  (B).  Nuclei  (N),  cell  membranes  (CM)  are  visible  in  the  tadpole 
images,  and  epidermal  cells  (EC)  can  be  seen  in  the  skin  images.  Images  were 
acquired  at  800  nm  wavelength  at  4  frames  per  second. 


with  a  confocal  gate  of  —30  pm.  While  this  confocal  section  thickness  is  nearly  6  times  that  of 
standard  histology,  a  high  transverse  resolution  of  <  2  pm  can  still  be  maintained.  Together  the 
combined  gating  effects  are  sufficient  for  high-resolution  in  vivo  imaging  of  cellular  features  in 
various  tissues,  demonstrated  in  Xenopus  laevis  tadpole  and  in  human  skin.  In  both  Xenopus 
tadpole  and  in  human  skin,  cellular  features  were  clearly  visualized. 
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4.1  OCM  Imaging  of  Pathology 


Pilot  ex  vivo  imaging  studies  were  conducted  in  the  clinical  pathology  laboratory  to  assess  the 
feasibility  for  cellular  resolution  imaging  in  endoscopy.  These  studies  were  performed  in 
collaboration  with  Dr.  James  Connolly,  M.D.,  Chief  of  Anatomic  Pathology,  from  the  Beth  Israel 
Deaconess  Medical  Center  in  Boston.  A  total  of  75  samples  from  39  patients  were  imaged. 
Normal  specimen  subtypes  included  squamous  esophagus  (10),  stomach  (3),  small  intestine  (5), 
colon  (11),  and  pancreas  (2).  Pathology  specimen  subtypes  included  columnar-lined  esophagus 
(11),  esophageal  adenocarcinoma  (1),  celiac  disease  (1),  inflammatory  bowel  disease  (3),  acute 
inflammation  (2),  chronic  colitis  (2),  melanosis  coli  (1),  tubular  adenoma  (11),  hyperplastic 
polyp  (5),  colorectal  adenocarcinoma  (4),  cholecystitis  (2),  and  chronic  pancreatitis  (1).  Figures 
4.4  and  31  present  sample  data  from  this  study.  Figure  4.4  a  illustrates  the  characteristic  pattern 
of  squamous  cells  in  the  esophagus.  Corresponding  histology  is  shown  in  Figure  4.4  b.  In  the 
OCM  images  as  well  as  histology,  cell  nuclei  can  be  clearly  differentiated  from  the  surrounding 
cytoplasm  and  individual  membranes  delineate  cell  boundaries.  Figures  4.4  c  and  4.4  d  show  an 
example  of  Barrett’s  esophagus.  Barrett’s  esophagus  is  a  condition  in  which  chronic 
gastrointestinal  reflux  leads  to  a  metaplastic  change  in  the  esophageal  mucosa  from  the  normal 
squamous  pattern  to  a  columnar  pattern  with  similar  features  to  colonic  mucosa.  The  presence  of 
Barrett’s  metaplasia  is  a  predisposing  risk  factor  for  the  development  of  dysplasia  and 
adenocarcinoma  of  the  esophagus.  The  hallmark  histopathologic  feature  of  Barrett’s  esophagus  is 
the  barrel-shaped  goblet  cell.  OCM  identifies  the  glandular  architecture  of  the  Barrett’s  mucosa 
as  well  as  the  presence  of  goblet  cells  in  the  columnar  epithelium. 


Figure  4.4.  OCM  images  and  histology  of  normal  esophagus  (A,B)  and  Barrett’s 
esophagus  (C,D)  in  vitro.  Normal  squamous  mucosa  has  squamous  cells  with 
centrally-located,  highly  scattering  nuclei  (n).  Barrett’s  epithelium  has 
intestinalized  glands  with  barrel-shaped  goblet  cells  (gc). 
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Figure  4.5.  OCM  and  histology  of  normal  (A,B)  and  dysplastic  (C,D)  colon  in 
vitro.  Normal  colonic  mucosa  has  round  crypts  with  goblet  cells  (gc)  and  basally 
situated  nuclei.  Adenomatous  dysplastic  crypts  have  more  eccentricity  and  cigar¬ 
shaped  nuclei. 

Figure  4.5  shows  a  comparison  of  OCM  images  of  normal  and  dysplastic  colonic  mucosa.  The 
normal  colonic  mucosa  shown  in  the  OCM  images  and  histology  of  Figures  4.5  a  and  4.5  b, 
respectively,  illustrates  a  regular  pattern  of  round  crypts  with  numerous  goblet  cells  and  nuclei 
restricted  to  the  basal  aspect  of  the  columnar  epithelium.  In  contrast,  Figures  4.5  c  and  4.5  d 
present  images  and  histology  from  a  tubular  adenoma  with  low  grade  dysplasia.  Glands  in  the 
adenoma  are  larger  and  exhibit  significant  eccentricity  compared  to  the  small  round  crypts  of  the 
normal  mucosa.  In  addition,  the  adenomatous  glands  show  the  presence  of  cigar-shaped  nuclei 
extending  beyond  the  basal  third  of  the  columnar  epithelium.  OCM  images  correlate  well  with 
histology  and  demonstrate  the  ability  to  identify  key  histologic  features  of  normal  and  pathologic 
tissues. 


4.3  Line  Scan  Full-field  OCM 

In  addition  to  continuing  research  on  OCM  using  standard  methods,  we  are  exploring  new 
techniques  which  can  enable  in  vivo  cellular  level  imaging.  OCM  is  typically  implemented  by 
scanning  a  tightly  focused  beam  in  a  raster  pattern  across  the  sample  to  generate  an  en  face 
image  at  constant  depth  [122,  124].  In  vivo  OCM  imaging  on  biological  samples  usually  requires 
the  imaging  speed  larger  than  several  frames  per  second  to  minimize  the  motion  artifacts. 
Therefore,  the  typical  raster  scanning  design  utilizes  fast  scanning  (-kHz)  in  one  dimension  (X) 
and  slow  scanning  (-Hz)  in  the  other  dimension  (Y).  A  significant  research  effort,  including  the 
work  described  in  the  above  section,  has  been  invested  on  developing  miniaturized  fast  scanning 
devices  for  XY  scanning  [125-127],  but  it  is  challenging  to  obtain  a  precise  beam  steering  and 
small  device  sizes.  Alternative  full-field  illumination  and  detection  approaches  have  also  been 
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developed  with  Linnik  type  interference  microscopes  which  eliminate  the  need  of  beam 
scanning,  and  therefore  enable  rapid  and  high  quality  image  acquisition  [128-134],  However,  the 
sensitivity  of  the  full-field  imaging  is  lower  compared  with  point  scanning  [135], 

We  developed  a  novel,  line  illumination  and  detection  scheme  to  perform  high-resolution  OCM. 
Line  scanning  removes  the  need  for  fast  raster  scanning  on  one  axis  (X),  therefore  allowing  the 
use  of  only  a  single  axis  slow  scanner  (Y)  in  the  sample  arm.  The  line  scan  OCM  system  is  based 
on  a  Linnik  type  interference  microscope  with  two  identical  microscope  objectives  placed  in 
both  the  sample  and  reference  arms,  as  shown  in  Figure  4.10  a.  A  compact,  commercially- 
available  broadband  Ti:Sappire  laser  was  used  to  generate  a  broadband  spectrum  around  800  nm. 
Line  illumination  was  achieved  by  using  the  combination  of  a  spherical  focusing  lens  and  a 
plano-concave  cylindrical  lens.  The  light  from  the  sample  and  reference  arms  are  recombined 
with  the  beam  splitter  and  imaged  onto  the  CCD  line  scan  camera.  The  signal  was  modulated  by 
moving  the  reference  mirror  with  a  piezoelectric  transducer  (PZT)  and  a  conventional  four-phase 
demodulation  was  employed  to  extract  the  interference  component  [129,  130].  En  face  images 
were  generated  by  transversely  scanning  the  sample  along  the  direction  orthogonal  to  the 
illumination  line  with  a  precision  translation  stage,  and  a  series  of  en  face  images  at  different 
depths  were  acquired  by  translating  the  sample  along  the  axial  direction. 


LSCCD 


Figure  4.10.  A.  Schematic  of  line  scan  OCM  system.  B.  Point  spread  function 
using  Ti:Sappire  laser  light  source.  C.  Image  of  resolution  chart.  Line  scan  is  a 
new  approach  which  promises  higher  speeds  and  reduced  scattering  compared 
with  full  field  OCT  techniques. 

In  the  preliminary  study,  the  line  scan  camera  was  operated  at  8  kHz  line  scan  acquisition  rate, 
and  the  PZT  was  sinusoidally  actuated  at  2  kHz  and  synchronized  by  the  camera.  10  lines  from 
each  of  the  four-quadrant  phases  were  accumulated  to  increase  the  signal-to-noise  ratio  and 
resulted  in  an  imaging  speed  of  200  transverse  lines/second.  Each  image  consisted  of  400 
transverse  lines  and  was  acquired  in  2  seconds.  The  field  of  view  (FOV)  on  the  camera  is  400 
pixels,  therefore  the  image  was  composed  of  400  by  400  pixels  with  the  size  of  400  pm  x  400 
pm.  Figure  4.10  b  shows  the  axial  resolution  of  the  system  measured  by  translating  the  mirror 
along  the  sample  arm,  and  a  ~3  pm  coherent  gate  is  achieved.  Figure  4.10  c  demonstrates  the 
coherence  gated  image  of  a  USAF  target  on  the  sample  arm.  The  smallest  Group  7,  Element  6 
bars  were  able  to  be  resolved,  which  shows  a  ~2  pm  transverse  resolution.  The  system  sensitivity 
was  measured  to  be  83  dB,  which  is  consistent  with  previous  publications  on  full-field  imaging 
[135]. 
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Line  scan  is  interesting  because  it  has  advantages  for  endoscopic  applications.  Although, 
transmission  of  the  line  information  through  the  endoscope  would  require  a  fiber  optic  ribbon,  a 
single  axis  slow  scanning  device  is  easier  to  implement  than  a  two  axis  device.  Line  illumination 
can  reduce  the  amount  of  incoherent  scattered  light  and  pixel  cross-talk  compared  to  full-field 
illumination,  therefore  would  improve  the  detection  sensitivity.  Furthermore,  the  use  of  line 
illumination  and  line  scan  cameras  promises  to  enable  reduced  sensitivity  to  sample  motion 
compared  to  full-field  detection  because  the  phase  sensitive  infonnation  is  acquired  more  rapidly 
before  phase  averaging  effects  can  occur. 

5.  Functional  Brain  Imaging  with  OCT 

Brain  neural  response  to  cognitive  and  sensory  stimuli  is  an  active  area  of  research.  Several 
imaging  techniques  have  been  applied  to  noninvasively  probe  neural  response  to  external  stimuli. 
Optical  imaging  is  uniquely  sensitive  to  neuronal  and  vascular  responses  to  brain  activation,  both 
invasively  [93,  136]  and  non-invasively  [94],  and  can  contribute  to  understanding  the  neuro¬ 
vascular  relationship  [137].  Currently  most  of  the  optical  imaging  presents  the  two-dimensional 
en  face  maps  of  brain  activation  [138],  Depth-resolved  images  can  be  achieved  by  optical 
tomography,  but  the  depth  resolution  is  limited  to  -100  pm  [139].  OCT  can  provide  depth- 
resolved,  micron-resolution  images  in  vivo  and  in  real  time.  OCT  can  also  provide  three- 
dimensional  volumetric  imaging  with  high  spatial  resolution. 

Several  studies  exist  to  date  which  use  OCT  to  study  functional  activation  in  neuronal  tissues. 
Maheswari,  et  al.  showed  depth  resolved  stimulus  specific  profiles  of  slow  processes  during 
functional  activation  in  the  cat  visual  cortex  [99].  They  attributed  these  signals  to  variation  in 
scattering  due  to  localized  structural  changes  such  as  capillary  dilation  and  cell  swelling. 
Lazebnik,  et  al.  subsequently  perfonned  a  simple  measurement  of  scattering  changes 
corresponding  to  fast  and  slow  signals  triggered  by  action  potential  propagation  in  the  sea  slug 
abdominal  ganglion  [97].  Satomura  et  al.  observed  the  delayed  swelling  of  the  cortical  surface  in 
the  somatosensory  cortex  following  the  electrical  stimulation  of  the  rat  hind  paw  [140],  and  Seki, 
et  al.  measured  the  cross-sectional  profiles  of  blood  flow  velocity  in  the  rat  pial  microvessels  and 
their  temporal  changes  [141].  Ooi,  et  al.  examined  the  neural  plasticity  in  response  to  peripheral 
neuropathic  pain  in  mice  in  vivo  using  OCT  [142].  Other  studies  have  also  investigated  precise 
measurement  of  nerve  axon  displacement  using  low  coherence  interferometry  methods  [143]. 

Our  collaborative  group  is  investigating  OCT  to  measure  subsurface  scattering  changes  due  to 
functional  activation  in  the  rat  somatosensory  cortex.  Simultaneous  OCT  and  video  microscopy 
imaging  were  performed  and  the  system  diagram  is  shown  in  Figure  5.1.  The  OCT  system 
operates  at  1060  nm  center  wavelength  using  an  Nd:Glass  laser  generating  -100  fs  pulses.  The 
laser  output  is  coupled  into  a  standard  HI- 1060  single  mode  fiber  and  the  nonlinear  effects  in  the 
single  mode  optical  fiber  broaden  the  optical  spectrum  to  -40  nm,  which  results  in  an  axial 
resolution  of  -18  pm  in  air  (corresponding  to  13  pm  in  tissue).  In  the  reference  arm,  a  delay 
scanning  at  -1150  Hz  is  performed  and  images  of  -380  transverse  pixels  are  acquired  at  3 
frames  per  second.  The  sample  beam  is  focused  to  a  spot  size  of  -40  pm  and  positioned  on  the 
rat  cortex  using  an  XY  scanning  pair  of  galvanometers.  The  system  allows  simultaneous  imaging 
with  the  video  microscopy  system.  The  measured  system  sensitivity  from  a  single  high  reflector 
was  >95  dB  with  -14  mW  illumination  on  the  sample. 
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The  video  microscopy  technique  illuminates  the  cortex  at  570  nm  and  images  scattered  light  onto 
a  two-dimensional  CCD  camera.  Changes  in  cortical  reflectivity  are  measured  in  response  to 
functional  activation.  The  reflectivity  changes  at  this  wavelength  are  primarily  due  to  change  in 
the  total  hemoglobin  concentration.  The  video  microscopy  method  is  used  in  this  experiment  to 
localize  and  characterize  the  functional  activation  region  for  placement  of  the  OCT  beam. 

Forepaw  stimulation  is  perfonned  using  20  second  stimulation  blocks.  Each  block  consists  of  1 
second  of  pre-stimulus  period  followed  by  4  seconds  of  stimulation  with  ~1.8  mA  pulses  at  3  Hz. 
A  15  second  post-stimulus  period  is  then  provided  to  allow  full  recovery  of  the  excitable  tissues 
to  baseline.  The  stimulus  block  is  repeated  60  times  during  data  acquisition  over  a  20  minute 
period.  Block  averaged  signals  are  then  computed  to  reduce  the  effects  of  physiologic  noise  in 
the  measurements.  For  data  processing,  digitized  OCT  images  are  converted  from  log  to  linear. 


Figure  5.1.  Schematic  diagram  of  simultaneous  OCT  and  video  microscopy 
system  for  measuring  functional  activation  in  the  rat  cortex.  The  forepaw 
stimulation  protocol  used  in  the  experiments  is  also  shown  here. 


All  experiments  were  conducted  according  to  protocols  approved  by  the  Committee  on  Animal 
Care  at  the  Massachusetts  Institute  of  Technology  and  the  Animal  Care  Committee  at  the 
Massachusetts  General  Hospital.  Rats  were  anesthetized  with  urethane  and  immobilized  in  a 
sterotaxis  before  beginning  the  cranial  window  preparation.  An  area  of  skull  overlying  the 
primary  somatosensory  cortex  on  the  contralateral  side  of  the  stimulated  forepaw  was  thinned 
with  a  dental  burr  until  transparent.  Thickness  was  measured  with  OCT  to  be  100-200  pm. 

Figure  5.2  shows  representative  OCT  imaging  of  rat  brain  structures.  OCT  cross-sectional 
images  are  exactly  co-registered  to  the  en  face  images  proved  by  video  microscopy.  The 
structures  in  rat  brain  including  the  skull  (S),  dura  (D),  cortex  (C),  and  subsurface  vessels  (V)  are 
discernable.  The  skull  appears  patchy  with  a  highly  scattering  layer  thought  to  be  the  dura 
separating  it  from  the  subarachnoid  space  and  the  cortex.  Large  surface  cortical  vessels  can  be 
seen  clearly. 
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Figure  5.2.  Left:  en-face  CCD  view  of  rat  brain  and  vascular  pattern  through 
cranial  window.  Red  arrow  denotes  the  OCT  imaging  plane.  Middle:  cross- 
sectional  OCT  image  of  rat  brain  showing  detailed  structures  including  skull  (S), 
dura  (D),  cortex  (C),  and  subsurface  vessels  (V).  Right:  schematics  of  brain 
anatomy. 

The  OCT  functional  signal  is  shown  in  Figure  5.3.  The  signal  is  computed  by  dividing  each  time 
point  in  the  stimulus  and  post-stimulus  periods  by  the  mean  of  the  pre-stimulus  period.  The 
signal  is  then  averaged  over  all  60  blocks  in  the  data  collection  period.  The  time  window  chosen 
for  plotting  Figure  5.3C  corresponds  to  mean  of  the  stimulus  period.  Several  local  regions  of 
interest  are  chosen  from  the  images  for  time  course  analysis.  These  regions  of  interest’s  appear 
as  boxes  of  different  colors  in  the  images  and  the  corresponding  time  courses  are  plotted  in 
Figure  3D.  The  time  courses  stay  at  baseline  during  the  pre-stimulus  period,  begin  rising  ~1 
second  after  the  onset  of  stimulus,  and  continue  to  rise  until  1-2  seconds  after  the  cessation  of  the 
stimulus.  The  time  course  returns  to  baseline  after  reaching  its  maximum.  The  response  time 
course  varies  in  amplitude  from  region  to  region.  For  the  most  part,  the  trend  of  the  response  is 
uniform  with  slight  variation  in  time  to  peak  response.  The  time  course  of  the  response  is  similar 
to  the  response  observed  with  the  video  microscopy  suggesting  that  the  OCT  signal  change  may 
arise  from  changes  in  total  hemoglobin  concentration. 
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Figure  5.3.  Representative  OCT  imaging  results  of  functional  activation  in  the  rat 
brain.  (A)  OCT  functional  activation  image.  (B)  Functional  signal  time  course  in 
region  of  interests  denoted  by  different  colors. 
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We  demonstrated  that  OCT  can  image  the  neuro-vascular  response  to  somatosensory  stimulation 
in  the  rat.  The  dominant  signal  grows  over  the  4  seconds  of  the  fore-paw  electrical  stimulus  and 
then  requires  ~5  seconds  to  recover  to  the  baseline.  The  time  course  of  the  OCT  response 
correlates  well  with  video  microscopy  suggesting  that  this  OCT  signal  may  reflect  changes  in  red 
blood  cell  density,  although  neuronal  swelling  can  not  be  ruled  out  at  this  point.  Localized 
regions  of  strong  signal  increases  where  identified  that  are  characteristic  of  blood  vessels.  In  the 
parenchyma,  smaller  signal  increases  and  decreases  were  observed.  This  is  consistent  with 
previous  observations  of  increases  and  decreases  in  red  blood  cell  density  in  the  capillaries 
during  brain  activation.  Further  work  is  necessary  to  confirm  the  origin  of  this  functional  signal. 
These  results  suggest  that  OCT  could  have  an  important  role  in  future  studies  of  the  functional 
neurovascular  response.  In  a  more  general  context,  the  integration  of  structural  and  functional 
imaging  in  ophthalmology  and  brain  sciences  promises  to  enhance  fundamental  understanding  of 
diseases  pathogenesis  as  well  as  develop  more  sensitive  markers  for  early  diagnosis  and 
assessment  of  disease  progression. 
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